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Lighthill's  theory  of  jet  noise,  as  extended  and  developed  by 
Ribner  (self  and  shear  noise),  has  successfully  described  many  featiires  of 
the  jet  noise  oirtside  the  'refraction  valley'.  However,  attesnpts  to  measvire 
the  self  and  shear  noise  soiirce  terms  directly  by  means  of  a cross- correlation 
technique  have  been  only  partially  successful.  The  major  difficulty  has 
been  suspected  as  spurious  'probe  noise'  generated  by  turbulence  - hot  wire 
interaction.  Thus,  to  avoid  this  problem,  the  traditional  hot  wire  anemometer 
has  been  replaced  in  the  present  investigation  by  a non- intrusive  device:  a 
Laser  Doppler  Velocimeter.  Substantial  modifications  were  smde  to  meet  the 
constraints  iii5>osed  by  the  correlation  experiment;  a major  feature  was  provision 
to  meas\ire  the  ccoponent  of  tvirbtilent  velocity  in  the  observer  direction  jc. 

Cross- correlations  and  cross- spectral  densities  of  the  jet  noise  at 
40"  to  the  jet  axis  and  the  instantaneous  turbulent  jet  flow  Ux  (^Ux/&t2  ~ 
shear  noise  source  term)  or  Ux^  (S^ux^/St^  ~ self  noise  source  term)  were 
measvired  at  various  source  positions  in  the  jet.  Source  distributions  were 
inferred  therefrom  over  slices  of  jet  normal  to  the  jet  axis.  They  were 
found  to  be  strongly  pear-shaped,  rather  than  axi symmetric,  the  small  end  of 
the  'pear'  pointiiig  toward  the  observer.  (This  is  not,  however,  inconpatible 
with  the  axlsymmetry  of  far  field  sound  intensity.) 

Self  and  shear  noise  spectra  have  been  constructed  from  the  measured 
cross- spectral  densities  by  a method  consistent  with  the  postulated  self/shear 
aolse  foimalism.  The  two  spectra  exhibit  conparable  amplitudes  and  virtually 
identical  shapes,  but  are  displaced  substantially  in  frequency:  all  this  is 
predicted  by  the  theory.  Self  and  shear  noise  spectra  extracted  from  far  field 
jet  noise  intensities  via  an  algorithm  of  Nossier  and  Ribner  exhibit  the  same 
behaviour.  On  the  whole  both  sets  of  spectra,  although  derived  from  vastly 
different  experimental  procedures,  are  conpatible. 
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1.  INTRODUCTION 


1.1  HlstorlceLl  Backeroiaid 

Over  a quarter  centviry  ago,  Sir  James  Lighthill  (Refs.  1,  2) 
proposed  a theory  of  sound  generation  by  turbulent  flows.  The  governing 
wave  equation  is  a direct  consequence  of  the  eqmtions  of  conservation  of 
mass  and  momentum.  The  sound  sources  are  described  as  an  equivalent  distri- 
bution of  quadrupoles  in  a qvilescent  medivmi  whose  strength  is  proportional 
to  time  varying  virtual  stresses  in  the  fluid.  For  a subsonic  Jet  at  ambient 
tenperature  the  inertial  Reynolds  stresses  pVj_Vi  are  the  dominant  sources. 
Scaling  considerations  led  Lighthill  to  deduce  the  u6d2  acoustic  power  law 
which  is  strongly  sipported  by  measxirements  from  subsonic  Jet  flows  (Ref.  3)- 

In  order  to  be  able  to  predict  the  relevant  properties  of  Jet  noise 
(the  mean  square  pressvtre  fot  exanple)  two  point  correlations  of  pv^Vi , or 
equlVEilent  soiirce  terms,  must  be  known  or  modelled.  Proudman  (Ref.  4;  pointed 
out  that  the  distribution  of  sources  with  strength  pvx^  (a  single 

term)  is  eqTiivalent  to  the  qxiadripole  distribution  of  strength  pvj^vj  (six 
terms) , if  the  radiation  to  an  observer  at  x far  away  from  the  source  region 
is  considered  (vx  is  the  velocity  in  the  direction  of  the  observer) . He 
applied  this  formalism  to  evaluate  the  noise  from  decaying  isotropic  turbulence 
described  by  a statistical  model. 


Lilley  (Ref.  5)  applied  Li^thill' s theory  in  a quantitative  fashion 
to  the  noise  generation  of  air  Jets.  A particvilar  featxire  of  his  approach  was 
a split  into  'shear  noise'  and  'self  noise'.  The  shear  noise  sources  were 
modelled  (Ref.  2)  as  proportional  t9  the  first  time  derivative  of  the  total 
pressure  multiplied  by  a mean  rate  of  strain.  Such  sources  were  believed  to 
be  dominant  in  regions  of  high  shear,  as  is  found  in  the  mixing  region  of  a 
Jet.  (Rlbner,  Ref.  6, plater  argued  that  the  relative  spatial  scales  of 
turbulence  and  mean  shear  Eire  ccanparable  in  such  a region  and  thus  the  shear 
noise  sources  aure  not  necessarily  domineuit.)  Additional  noise  - self  noise  - 
is  generated  by  the  turbulent  flow  regardless  of  the  presence  of  mean  shear. 
Lilley  ccoputed  hypothetical  spectra  and  directivity  patterns  based  on  assumed 
two  point  correlations  of  the  self  and  shear  noise  soxirces.  He  independently 
derived  the  x°  and  x""^  laws,  first  obtained  by  Ribner  (Ref.  7),  which  describe 
the  axial  source  strength  distribution  in  the  mixing  region  and  the  fully 
developed  Jet  respectively.  Powell  (Ref.  8)  also  developed  these  scaling 
laws  from  tiirbulence  similarity  considerations,  as  well  as  a corresponding 
Jet  noise  spectral  shape  (f2  at  low  frequencies,  f"^  at  high  frequencies). 


A different  self  and  shear  noise  model  was  proposed  by  Ribner 
(Refs.  9>  10)  in  which  the  momentum  flux  pvx^  of  the  Proudman  formalism  is 
decon5)08ed  into  two  contributing  terms  2pUxibc  [Ux(^n3{/^t2)  ~ shear  noise] 
and  pUx^[02ux^/St^)  ~ self  noise].  He  was  able  to  predict  the  qualitative 
behaviour  of  the  self  and  shear  noise  both  as  to  directivity  and  relative 
pitch.  Source  convection  was  included  in  the  analysis,  and  good  agreement 
with  measurement  was  found  for  field  points  outside  the  valley  of  relative 
quiet  that  exists  near  the  Jet  axis. 


The  earlier  work  (e.g.,  Lighthill,  Lilley,  loc  clt)  Included 
source  convection  (at  a fraction  Me  of  the  speed  of  sound)  but  neglected 
small  variations  in  the  travel  time  (the  time  delay  between  emission  and 
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reception)  across  a coherent  source  region;  it  predicted  a (l  - M^cose)”^ 
aii5)lifying  'convection  factor'  which  can  be  singular  for  Me  > !•  Ffowes 
Williams  (Ref.  U)  and  Ribner  (Ref.  6)  showed  by  different  methods  that  the 
singxxlarity  is  removed  if  the  correct  time  delay  is  used.  According  to  this 
factor,  sovirce  convection  at  subsonic  speeds  (Me  < 1)  strongly  biases  the  noise 
towards  the  downstream  direction  (0  -+0).  However,  refraction  bends  the  sound 
rays  towards  the  sides  to  create  the  peciiliar  'valley  of  quiet'  along  the  jet 
axis;  this  gives  the  overall  directivity  a heart-shaped  ^pearance. 

The  role  of  refraction  by  the  jet  velocity  and  teii5)erature  gradients 
was  recognized  independently  in  the  United  Kingdom  (Ref.  12)  and  in  Canada. 
Ribner 's  grovq),  in  a major  series  of  e^eriments  (Refs.  13,  1^),  showed  that 
the  soiuid  of  a pure  tone  point  sovirce  placed  in  a jet  flow  has  a valley  of 
quiet  that  at  each  frequency  closely  matches  that  of  narrow  band  filtered  jet 
noise  centred  at  the  same  frequency,  ffiince  the  pure  tone  source  is  omnidirec- 
tional when  the  jet  is  turned  off,  this  confirms  that  the  valley  of  quiet  in 
the  jet  noise  must  be  due  to  refraction  of  sovmd  by  the  mean  flow  velocity- 
ten?)erature  gradients.  Later  Schubert  (Refs.  15,  l6)  theoretically  predicted 
similar  directivities.  The  noise  in  the  'valley  of  quiet'  is  low  pitched;  the 
pitch  increases  steadily  as  the  observer  moves  away  from  the  jet  axis;  this  ' 
is  just  the  reverse  of  what  would  be  expected  from  simple  Doppler  shift 
considerations.  However,  high  frequencies  suffer  more  refraction  than  low 
frequencies  and  are  'bent'  further  away  from  the  jet  axis;  this  in  conjunction 
with  the  directional  and  pitch  properties  of  Ribner 's  self  and  shear  noise 
model  explains  the  'reverse  Doppler  shift' . 

MacGregor  et  al  (Ref  17)  demonstrated  a viable  semi-en^iirical  theory 
based  on  Ribner 's  self  and  shear  noise  model  (which  neglects  refraction)  in- 
corporating measured  or  calculated  (Refs.  15,  l6)  refraction  corrections.  The 
en^iiricism  of  the  model  correlations  is  a substitute  for  a detailed  spatial 
and  ten^ioral  knowledge  of  the  relevant  turbulence  parameters,  and  the  a poster- 
iori refraction  correction  aJ-lows  for  reasonable  estimates  in  the  refraction 
zone . 

Refraction  is  automatically  accounted  for  if  a convected  wave 
equation  with  variable  mean  flow  is  considered.  An  equation  of  this  kind 
was  derived  by  Phillips  and  applied  to  predict  the  sound  from  a turbulent 
supersonic  shear  flow  (Ref.  l8) . A similar  equation  has  been  used  in  the  study 
of  combustion  noise  (e.g..  Ref.  19) , although  with  different  dominant  source 
terms.  Lilley  (Ref.  20)  went  furtiier  and  formulated  a third  order  generalized 
convected  wave  equation  which  contrasts  with  the  more  traditional  second  order 
equations  of  the  previous  models.  The  Lilley  equation  sepsirates  generation 
effects  (sources)  from  propagation  effects  more  precisely  than  the  Lighthill 
formalism;  thus, despite  its  complexity,  it  has  come  into  widespread  use  (Refs. 
21,  22). 

The  convected  wave  equation  can  only  be  solved  analytically  for 
idealized  mean  flow  distributions.  Mani's  (Ref.  23)  model  of  jet  noise 
generation  is  one  of  them.  It  use*  a form  equivalent  to  Lilley' s equation, 
and  describes  moving  quadrupoles  in  a cylindrical  jet  embedded  in  a fluid  at 
rest.  In  Mani's  model  the  shear  term,  which  in  Ribner 's  model  is  on  the  right 
hand  side  of  the  wave  equation  as  a source  term,  is  in  effect  on  the  left  hand 
side  as  a prepagation  term.  Ribner  points  this  out  in  Ref.  2k  and  demonstrates 
that  the  respective  theories,  although  they  differ  in  the  role  of  the  shear 
term,  show  comparably  good  agreement  with  experiment  outside  the  refraction 
valley. 
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Recently  recognition  of  the  existence  of  large-scale  coherent 


structures  in  jet  tixrbulence  has  led  to  analyses  (Refs.  24,  26)  which  suggest 
that  such  structures,  particularly  axi symmetric  ones,  can  be  efficient  soiuid 
radiators.  The  experimental  evidence,  however,  fails  to  support  this  in  the 
case  of  jet  noise.  Extensive  far-field  measurements  of  Maestrello  (Ref.  27) 
show  that  two  microphones  on  opposite  sides  of  a jet  have  very  poor  signal 
correlations  except  at  small  angles  to  the  axis;  thus  any  axisymmetric  com- 
ponents contribute  only  a gmal l fraction  to  the  total  jet  noise.  Furthermore, 
Ribner  (Ref.  28)  has  recently  shown  that  his  self  and  shear  noise  model  is 
capable  of  predicting  two  point  correlations  of  the  jet  noise,  in  good  agree- 
ment with  measurements  by  Maestrello,  without  invoking  large  scale  structures. 

An  altogether  different  approach,  introduced  circa  1958,  models  the 
soiurces  as  flxiid  dilatations  proportional  to  fluctuating  pressures  (pseudosound) 
in  and  near  the  jet  flow.  Such  a source  model  is  inplicit  in  the  work  of 
Meecham  and  Ford  (Ref.  29);  Corcos  and  Broadwell  (Ref.  30)  an<i  Ribner  (Ref.  31) 
formulated  the  model  explicitly.  Ribner  (Ref.  6)  developed  the  dilatatior 
approach  in  detail  and  illustrated  phenomena  associated  with  source  convection 
and  wave  cancellation  effects,  among  others,  all  of  which  can  also  be  applied 
to  Lighthill’s  formalism.  The  equivalence  of  the  dilatation  and  Lighthill 
models  has  been  demonstrated  by  Ribner  (Ref . 6)  and  acknowledged  by  Lighthill 
(Ref.  32). 

At  first  the  dilatation  theory  did  not  enjoy  widespread  attention, 
possibly  for  lack  of  detailed  experimental  information  about  the  pseudosoimd 
field  of  a turbulent  jet.  On  the  experimental  side  there  were  suspicions  that 
pressure  probes  locally  distort  the  pseudosound  pressure  they  were  intended 
to  measure.  Over  the  years  the  experimental  difficulties  have  been  overcame, 
and  the  dilatation  theory  continues  to  be  the  subject  of  experimental  and 
theoretical  investigations  (Refs.  33}  34,  35 > 36). 

1.2  Ribner ’s  Self  and  Shear  Noise  Model 

Ribner  (Refs.  9,  10)  developed  the  Proudman  formalism  by  noting 
that  the  instantaneous  velocity  in  the  direction  of  an  observer  at  x in  the 
far  field  is  the  sram  of  a mean  velocity  Ux  and  a turbrjlent  velocity  Hv(y,t)  . 
This  leads  to  two  families  of  sources:  one,  proportional  to  contains 

only  the  rmsteady  turbrolent  velocity.  Noise  generated  by  such  sources,  which 
are  found  in  all  turbulent  flows,  is  called  self  noise.  For  the  case  of 
isotropic  turbxilence  (assumed  in  his  analysis)  sources  of  the  type  2Ux(8^x/^^) 
can  be  shown  to  contribute  only  if  the  mean  flow  possesses  shear;  hence  the 
name  shear  noise  is  associated  with  this  second  source  term. 

Using  certain  mean  flow  and  turbrolence  models,  Ribner  was  able  to 
predict  the  qualitative  behaviour  of  the  self  and  shear  noise.  The  results 
suggest  that  the  basic  self  noise  is  omnidirectional,  whereas  the  shear  noise 
has  a dipole-like  directivity  with  lobes  pointing  along  the  jet  axis.  The 
shear  noise  spectrxim  is  predicted  to  peak  about  one  octave  below  the  self  noise 
spectrum;  the  maximimi  Intensities  are  comparable.  These  basic  patterns  are 
powerfully  altered  by  convection  and  refraction  effects. 

The  first  extensive  comparison  of  the  self  and  shear  noise  model 
with  experimental  data  was  performed  by  MacGregor  et  al  (Ref.  17).  More 
recently  Nossier  and  Ribner  (Ref.  37)  found  good  agreement  between  predicted 
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and  measured  narrow-band  filtered  jet  noise  directivity  patterns  and  demon- 
strated that  the  self  and  shear  noise  spectra  extracted  from  far  field  data 
of  several  investigators  conformed  to  the  predictions  of  the  theory.  The 
conjjari sons  were  restricted  to  points  outside  the  refraction  valley,  as  the 
basic  theory  does  not  predict  refraction. 

The  surprising  lack  of  spatial  coherence  of  the  far  field  jet  noise 
inferred  from  measured  two  point  correlations  of  jet  noise  (Maestrello,  Ref. 

27)  led  Ribner  (Ref.  28)  to  extend  his  self -shear  noise  pre^ctive  model  to 
deal  with  two  field  points  (microphones)  rather  than  one.  The  result  was  a 
quite  reasonable  quantitative  fit  of  the  normalized  theoretical  arid  measured 
cross- correlations.  The  agreement  is  especially  good  for  the  plane  at  90® 
to  the  jet  axis  \diich  should  be  insensitive  to  convection  effects  (the  model 
excluded  convection  effects  as  a concession  to  mathematical  expediency) . 

At  the  present  time  this  is  the  only  theory  that  has  been  demon- 
strated to  account  for  the  major  features  of  jet  noise  for  both  one  and  two 
microphones  outside  the  zone  of  refraction.  The  model  contains  one  empir- 
ical .constant;  rather  idealized  correlations  are  used  and  the  average  temporal 
scale  (or  equivalently  the  Strouhal  number)  is  chosen  for  best  average  fit; 
the  Stroiohal  number  is  roughly  conpatlble  with  measirrements  in  jet  turbulence. 

It  can  be  argued  that  Lighthill's  theory,  which  underlies  all  this, 
describes  only  an  equivalent  soiirce  distribution  and  not  the  'reail'  one,  since 
the  source  integral  has  been  modified  by  the  application  of  two  integral 
transformations.  While  the  value  of  the  integrals  may  well  be  equal,  the 
integrands,  which  are  proportional  to  the  effective  source  strength  distri- 
butions, may  not  be  one  to  one.  The  dilemma  is  analogous  to  the  one  foimd 
in  electrostatics,  where  dissimilar  distributions  of  monqpoles  and  dipoles 
can  resiilt  in  equivalent  electric  fields.  (A  monopole  distribution  of  local 
strength  equal  to  the  minus  divergence  of  a certain  dipole  distribution  will 
generate  the  same  field.) 

Even  though  Ribner' s self  and  shear  noise  formalism  appears  to  be  a 
valid  description  of  noise  generation  by  turbvilent  flows,  there  has  been  a 
lack  of  conclusive  evidence  such  as  furnished  by  direct  measurement  of  the 
supposed  so\irces  leading  to  correct  prediction  of  the  far  field  noise.  The 
first  serious  attenpt  to  measure  the  self  and  shear  noise  was  made  by  Chu 
(Ref.  38);  he  measiired  pertinent  two  point  velocity  correlations  in  a jet 
flow.  The  fourth  time  derivative  of  these  correlations  can,  in  principle,  be 
used  to  determine  the  properties  of  the  radiated  noise,  and  perhaps  lead  to 
the  validation  of  the  self  and  shear  noise  model.  It  was  found  that  multiple 
differentiation  of  time  domain  data  is  invariably  acconpanied  by  a large  error. 
Later  Lee  and  Ribner  (Ref.  39)  used  an  idea  of  Siddon  to  estimate  the 
combined  self  and  shear  noise  source  strength  distribution  as  well  as  the  far 
field  spectr\mi  from  two  point  cross-correlations  between  measurements  of 
Yx^  = (Ujc  iix)^  and  radiated  sound.  Only  two  differentiations  are  needed; 
however,  there  is  a trade-off  between  the  relative  ease  of  data  analysis  and 
poor  signal  to  noise  ratio,  as  the  cross- correlations  by  natiire  possess  very 
■small  ajiplitudes.  Experimental  difficulties  finally  led  Lee  to  abandon  the 
correlation  in  favour  of  the  mathematically  eqxd. valent  cross- spectral  density. 
With  the  latter  approach  estimates  of  the  overall  sotirce  strength  distribution 
and  the  radiated  spectrum  were  obtained. 
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Seiner  and  Reethof  (Ref.  4o)  applied  this  technique,  but  with  the 
hot  wire  signal  split  into  IJx  and  Ux  to  extract  the  apparent  self  and  shear 
noise  source  strength  distributions.  Their  measvireinents , supported  by 
earlier  unpiiblished  work  of  Morris  in  this  laboratory  (Ref.  4l) , appeared  | 

to  show  that  the  shear  noise  term  doninated  the  self  noise  by  at  least  10  dB, 
whereas  the  theory,  supported  by  much  indirect  evidence  from  far  field  jet 
noise  measurements,  sviggests  that  the  self  and  shear  noise  should  be  comparable 
in  intensity.  It  was  this  discrepancy  that  provided  much  of  the  motivation  I 

for  the  present  in-depth  investigation. 

i 

In  the  direct  correlations  of  the  jet  flow  and  jet  noise  discussed  I 

above,  the  instantaneous  jet  velocity  components  Ux  and  Ux  were  measured  with 
hot  wire  probes.  The  presence  of  such  solid  objects  - the  probe  support 
needles  - in  a turbxilent  flow  gives  rise  to  noise  generated  by  the  unsteady  j 

forces  that  act  thereon  (Refs.  42,  43)-  Even  though  the  very  weak  sound 
radiated  by  the  hot  wire  probe  is  biaried  in  the  overall  jet  noise,  the  small 
signal  may  be  well  correlated  with  the  flow  that  impinges  on  the  probe  and 
is  measured  by  it.  This  implies  that  the  measured  correlations  of  the  probe  j 

signal  and  the  far  field  microphone  signal  could  be  significantly  ' contaminated'  { 

by  a spurious  probe  noise  contribution.  This  second  possibility  further 
helped  motivate  the  present  experimental  investigation.  (The  theoretical 
expectations  for  probe-flow  interaction  are  discussed  in  detail  in  Appendix  B.)  j 

1.3  Present  Investigation  ! 

The  noise  generated  by  the  probe-flow  interaction  will  be  eliminated 
when  the  hot  wire  anemometer  is  replaced  by  a remote  sensing  scheme  such  as 
the  Laser  Doppler  Velocimeter  (iDV) . Since  the  demonstration  of  the  Laser 
Doppler  method  by  Yeh  and  Cummins  in  1964  (Ref.  44)  , the  technique  has  been 
refined,  and  reliable  Laser  Doppler  systems  have  been  available  'off  the  shelf 
for  several  years. 

Knott  et  al  (Ref.  45)  have  performed  parametric  LDV  studies  of  a 
high  speed,  high  temperatiire , turbulent  jet,  and  eilso  demonstrated  that  the 
LDV  can  be  used  in  cross- correlation  measurements  of  the  jet  flow  and  jet  noise. 

The  cross- correlation  measiirement  was,  however,  restricted  to  a single  pair  of 
field  points  in  the  jet  flow  and  the  far  field. 

The  present  investigation  is  an  extension  of  the  work  of  Lee  and  j 

Ribner  (Ref.  39)  and  others  (Refs.  40,  4l)  in  the  field  of  jet  noise  diag-  ^ 

nostics  via  cross- correlations . A major  difference  is  the  replacement  of  the 
hot  wire,  with  its  spiirious  noise,  by  an  LDV.  The  resesirch  program  was 
divided  into  three  tasks: 

(i)  adaptation  of  an  'off  the  shelf  LDV  system  to  the  special 
requirements  of  a cause  and  effect  experiment, 

(ii)  identification  of  possible  sources  of  error  that  may  affect 

the  estimates  of  the  correlation  between  jet  flow  and  jet  j 

noise,  as  well  as  validation  of  the  Laser  Doppler  Velocimeter 
system, 

(iii)  direct  correlation  of  the  postulated  self  and  shear  noise 
soiirce  terms  at  selected  field  points  in  the  jet  with  the 
reidlated  sound,  sind  ccB5)utation  of  overall  self  and  shear 
noise  spectra  based  on  the  correlation  measurements. 
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Assessment  of  the  coiq>atiblllty  of  the  predicted  fax-field  spectra  with 
measured  curves  would  test  whether  the  self  and  shear  noise  model  is  a self- 
conslstent  descrlptloi  of  sound  generated  by  turbulent  flows. 


2.  THEORY 

2.1  Fundamental  Equations 

The  central  equation  of  noise  generation  by  tuibulent  flow  can  be 
derived  from  the  equations  of  mass  and  momentm: 


= 0 


(2.2) 


p is  the  fluid  density,  P the  pressxire,  v^  the  fluid  velocity  in  the  i direction, 
and  Tij  the  viscous  stress  tensor.  (i,j)  range  from  1 to  3>  and  r^eated  indices 
are  sunned  over.  The  operation: 


+ v^I^) 


(2.3) 


combines  Eqs.  2.1  and  2.2  into  the  form: 


1 

Cq  at 


1 


Co^  at 


(p 


(2.4) 


after  sidding  Cq"^  a^/at^(P)  to  both  sides  (see  for  exanqple  Ref.  6).  The  left 
hand  side  of  Eq.  2.4  is  the  linear  wave  operator  applied  to  the  pressure  P(y,t) . 
In  a ccopresslble  stationary  medium,  small  disturbances  propagate  at  a speed 
Co  that  is  characteristic  of  the  medium  (Ref.  46).  Equation  2.4  therefore  can 
be  thought  of  as  describli^g  the  soiaid  field  of  a turbulent  flow  by  an  equivalent 
distribution  of  sound  sources  in  a quiescent  medixan  (Ref.  l)  with  no  flow  and 
uniform  sound  speed  Cr>. 


The  sinplicity  of  the  equation  is  deceiving,  as  it  allows  for  all  the 
subtleties  of  sovirce  motion,  refraction  of  sound  by  mean  flow  and  tenperature 
gradients,  and  other  phenomena.  Pertiubations  in  the  density  p for  example,  in 
conjunction  with  mean  velocity  gradients  are  responsible  for  refraction  of 
sound.  If  the  flow  is  anisentropic  (e.g.,  a heated  Jet),  the  term  c©"^ 

(P  - Cq^p)  will  be  nonzero.  This  causes  further  refraction  and  can  be  a signifi- 
cant sovurce  of  sound  (combustion  noise  for  example).  There  is  no  unique  descrip- 
tion of  sound  generation  by  a turbulent  flow,  as  the  conservation  equations  and 
their  solutions  may  be  manlpiU.ated  in  a variety  of  ways  (Refs.  l8,  20). 

In  Llgh thill's  treatment  of  a turbulent  subsonic  Jet  ^ ambient 
temperature  the  sovirce-like  terms  are  taken  to  be  dominated  by  S^PoVi’^j/SoTiSiyj  > 
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and  the  features  mentioned  above  are  sxjppressed.  The  principal  contributor 
to,  this  svH>ressicii  is  the  replacement  of  the  variable  local  density  p by 
the  constant  ambient  density  Pq.  Thus,  as  a first  approximation,  the  sound 
generated  by  such  a flow  is  described  by  the  differential  equation; 


^=0 


5t‘ 


^2 

o p V.v. 
O i j 


(2-5) 


The  pressure  at  any  field  point  x and  time  t is  confuted  by  the 
volume  integral  (Ref.  4?): 


= ^o  ^ feJJ 


k-  il 


<2-6) 

^ J 


The  source  region  is  assumed  to  extend  over  a volume  V that  encloses  the  region 
of  sensible  (i.e.  measurable)  ^jet  flow.  The  integral  sians  the  contribution  of 
each  elementary  sovu'ce  at in  V.  The  delta  function  determines  the  difference 
time  of  signal  emission  at  y and  signal  reception  by  an  observer  at 
is  the  source- receiver  separation. 


between  the 
X.  U - y| 


A mathematically  equivalent  description  is  given  by; 

1 rr  r - y|) 


dt'dy 


ffj- 

- 6(t'  - t + c"-^  X - y|)  - 

Jj  ^ 

L Ix-yl  J 

Application  of  the  divergence  theorem  (Gauss  theorem)  shows  that  the  first 
integral  does  not  contribute  to  P(x,t),  if  there  are  no  solid  boundaries:  the 
transform  into  an  integral  over  th¥  bounding  surface  of  V is  zero,  as  the 
integrand  decays  faster  than  r"^  for  large  r.  The  procedure  may  be  repeated 
for  the  remaining  integral: 


^2  6(t’  - t + c"^lx  - y]) 


I*  - yl 


dt'dy  (2.8) 


where 


^ <=o^k-  yl) 


I*  - y| 

1 r (■‘i  - yjXx.,  - y,) 


-1 


at' 


6(t'  - t + c^  1*  - yl) 


Sjj  - yjX^.i  - y.i)  a . . _-i 


k-zl 


IS  - y| 


- gfT  6(f  - 1 + =;  I*  - y|) 


— ■= 6(t’  - t + c 

I*- y| 


"li-y|>) 


(2.9) 
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and 


^2  5(f  - t + c"^|x  - y|) 

^ Ix-yl 


- - y|); 


|x|»|yl  (2.9A) 

"»  -7 


For  an  observer  outside  the  flow  and  sufficiently  far  away  so  that 
lx  - jr|  **»  [x|  » typical  wavelength,  the  only  significant  contribution  to  the 
acoxistic  pressxire  is  made  by 


p(4,t)  - p„ 


p v.v. 

^O  1 J 


6(t’ 


""  (2.1o7 


since  the  other  terms  decay  faster  than  12^1"^. 

Noting  that  xiVi|x|“^  = ^ • v|x|"l  = vx  the  conponent  of  the  velocity 
in  the  direction  of  the  observer  at  x,  and  integrating  by  parts  results  in  the 
so-called  Proudman  formalism: 


p(x,t)  = P(x,t)  - P^  ^ 

4t7C^  X 


:fj ^ ^ - yl^ 

I V (2.11) 


In  this  particular  description  only  the  velocity  component  in  the  direction  of 
the  observer  contributes  to  the  sovind  heard  by  the  observer.  Vx  can  be  measured 
with  considerably  less  effort  than  v^vj  (simultaneous  measurement  of  v^vj  is 
a Herculian  task),  and  thus  opens  the  way  for  experimental  verification  of  the 
theory  via  a correlation  technique. 


Equations  2.10  and  2.11  are  the  central  equations  governing  the 
generation  of  sound  by  unsteady  flows;  they  are  referred  to  as  the  Lighthill 
and  Proudman  forms,  respectively. 


2.2  Auto-Correlation  and  Power  Spectral  Density  of  the  Far  Field  Pressure 


The  instantaneovis  far  field  acoustic  pressure  (the  pressure  differ- 
ence from  ambient)  due  to  a turbulent  Jet  flow  fluctuates  above  and  below  zero 
in  a random  manner.  It  is  therefore  appropriate  to  characterize  the  noise  by 
measuring  relevant  statistical  properties  such  as  the  auto- correlation  function, 
and  the  power  spectral  density  of  the  pressure.  The  auto- correlation  of  the 
far  field  acoustic  pressure  is  defined  to  be  (Ref.  48) : 


R t,(xjt)  = p(x,t)  p(x,t  - t)  = 
pp  -r  ' -r  ' T -+«» 


substituting  for  p(x,t): 


p(x,t  - t)  dt  (2.12) 

- t + c"^  [x  - y I) 


p(x,t  - T)dt'dt(ty 
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(2.13) 


In  general  the  velocity  Vx  is  the  svm  of  a mean  velocity  Ux  ® ^x  and  a turbulent 
velocity  = Vx  - Vx*  The  source  term  expands  into 


(u. 


2U  u 

X X 


u ) 


(2.14) 


The  first  term  is  identically  zero,  hence; 

2 " 


f ^ f =0^  Is.  - Jjl)  S (^-^5) 

L X vJ.  dt  ^ 


^’^o  1^1  V 


For  steady  Jet  flow  u^j.,  and  p are  statistically  stationary  random 

variables  (Ref . 49) , so  that 


27rc 


T, — f ^x^^^  ^ \(y>4)  p(x,t  - T + c"^|x  - y|) 
X J * dr  ^ " o ^ 


o 1^1  V 


(^>t)  p(x,t  - 


I)  dy 


(2.16) 


The  far  field  acoustic  pressure  auto- correlation  is  the  sum  of  two 
volume  integrals.  At  t = 0,  the  integrals  describe  the  mean  square  pressxire. 
Ribner  (Ref.  10)  has  called  the  two  contributions  shear  and  self  noise  respec- 
tively. This  kind  of  division  was  first  suggested  by  Lilley  (Ref.  5)>  but 
with  different  definitions. 

Following  the  notion  of  Siddon  (Ref.  50),  one  may  think  of  these 
cross-correlations  as  causality  correlations , as  they  compare  the  st5)posed 
source  terms  [2Ut  ^/^^^(ujj),  d^/dt2(u3j2) ^ with  the  effect  they  produce.  The 
correlations  have  dimensions  of  a mean  square  pressure  per  unit  volume  and 
describe  the  sound  radiated  by  a correlation  volume  about  the  soizrce  point. 

A single  correlation  volimoe,  therefore,  contributes  a amai i amount  to  the  mean 
square  pressttre.  If  one  assianes  that  there  are  N such  correlation  volumes  of 
equal  strength  (AV*N  = V),  then  it  follows  that  the  normalized  peak  correlation 
coefficient  is  of  order  i/*»^  (see  Appendix  A for  a more  detedled  mathematical 
analysis) . 

The  Wiener- Khlnchin  theorem  states  that  the  auto- correlation  function 
and  the  power  spectral  density  are  Fourier  transform  pairs  (Ref.  51); 

00 

= J'  4>pp(x>w)  dw  (2.17) 
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The  power  spectral  density  cppp(x,w)  can  be  determined  from  Eq.  2.12  by  taking 
the  inverse  Fourier  transform; 

Substituting  for  F.pp(x,T): 

'Ppp^^’")  = Ux(y)  ^ u^(y,t)p(x,t  - T + c;^lx  - ^1)  6^“"^  dT<^ 


^T-Vr-riTrs  - T + c;^lx  - y|)  e^""^  dTdy  (2.19) 

8v  c |xlJJ  &T 


o '-t'  V 
After  integration  by  parts 

P, 


iu)c;^  Ix-y  1 


dy 


p 2 iuc"^  [x-y  I 

— u)  cp  p (x,y,a))  e ° dy 

X J u p 


kyrc 


(2.20) 


O W V 


where 


^u  "hf  ^TyTtrpTxrr^  e^“-'  dr 


(2.21) 


and 


9 2 (2.22) 

'^x  ^ 

The  power  spectral  density  of  the  far  field  sound  pressure  is  expressed 
by  Eqs.  2.20  to  2.22  as  a volume  integral  of  two  wei^ted  cross- spectral  density 
functions  9ujjp(x,y,w)  and  9uj^2p(x,y,u)) . The  double  differentiation  with  respect 

to  T in  the  time  domain  is  equivalent  to  multiplication  by  -u^  in  the  frequency 
domain  (Ref.  52),  a property  that  will  be  used  in  the  analysis  of  the  measured 
data. 


The  cross- correlation  fvmctions  Hujjp(x,y,T)  = Ux(y,t)  p(x,t  - t)  and 

Ru3j2p(5>y>-r)  = Ux2(^,t)  p(x  ,t  - t)  as  well  as  the  cross- spectral  densities 

9uj^(^>y»“)  ‘Pujj2p^4’y’‘^^  estimated  by  direct  measurement.  Both 

representations  contain  the  same  information,  and  it  is  often  a matter  of 
convenience  which  one  is  to  be  used.  For  each  position  x the  above  functions 
must  be  defined  over  all  values  of  y,  t,  or  w. 
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As  turbulence  levels  vary  slowly  with  position  in  the  jet,  the  volxnne 
integrals  2.l6  and  2. 2d  may  be  approximated  by  a suitably  weighted  sunoation 
over  a finite  nvi]!i>er  of  source  points.  The  welj^tlng  is  a measure  of  the 
effective  volume  over  \dilch  the  single  estimate  is  a reasonable  approximation 
to  the  integrand.  For  exanple,  the  ejqpression 


N 


R (x,t) 

PP  v’  • 


Situ, 


i7 


(2.23) 


estimates  the  auto- correlation  of  the  acovistic  pressure  from  measured  'causality 
correlations'  of  the  postulated  sources  (as  described  by  Ribner's  self  and  shear 
noise  formalism)  and  the  far  field  pressure. 


The  time  domain  estimates  are  very  susceptible  to  signal  to  noise 
problems  becavise  of  the  aman  aaplitudes  of  the  cross- correlations.  On  going 
over  to  the  frequency  domain,  the  transform 


/ N . ^o  ^ 2 , , ^"®o  l2S^JfnL„(l) 

^ - —rr,  A V?n^“  'Pu  ^n 

2t^o  n=l  ^ 


o 


hjic  ixl  up 

o '-7'  m=l  X ^ 


M , -1 1 I 

V 2 , ^ 

) U qp  2 (^>yTn»“^® 

A u p " 


(2.24) 


describes  the  power  spectrum  of  the  fax  field  pressure  'as  a simmation  of 
measured  cross- spectral  densities.  The  differentiations  with  respect  to  t 
of  the  time  domain  data  will  invariably  result  in  more  errors,  whereas  multi- 
plication by  -w^  should  only  alter  the  shape  of  the  cross-spectral  densities, 
which  are  likely  to  be  smooth  functions  of  frequency. 


The  measurement  of  the  'causality'  correlations  and  the  cross- spectral 
densities  is  one  of  the  major  objectives  of  this  investigation.  Should  the 
experimental  results  be  demonstrated  to  be  compatible  with  the  above  eqizations, 
then  the  self  and  sheau  noise  model  becomes  a more  credible  description  of  the 
noise  generation  of  tvirbulent  Jets. 

It  must  be  emphasized  that  the  measured  cross-correlation  and  cross- 
spectral  densities  are  not  estimates  of  a unique  source  distribution;  they  are 
but  estimates  of  one  of  a ntimber  of  mathematically  equivalent  source  patterns. 
The  overall  contribution  to  the  far  field  sound  pressure  from  such  eqvdvalent 
source  distributions  are  identical. 
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The  relative  Intensity  and  directivity  of  the  self  and  shear  noise 
have  been  predicted  theoretically  (Ref.  10)  \islng  certain  models  to  describe 
the  statistics  of  the  turbvdence.  The  self  noise  (aasisning  isotropic  tur- 
bvilence)  is  omnidirectional  and  the  directional  pattern  of  the  shear  noise  is 
dipole-like,  with  the  lobes  pointing  along  the  jet  axis.  The  basic  overall 
directivity  is  significantly  modified  by  source  convection  and  refraction 
(Refs.  13,  l4,  15,  16).  Refraction  of  soiind  by  a jet  flow  is  generally 
significant  only  for  angles  less  than  4o®  from  the  jet  axis.  Any  experiment 
designed  to  detect  a cause-effect  relationship  between  the  principal  source 
term  and  the  radiated  far  field  sound  pressure  should  avoid  measvirements  in 
that  region,  unless  it  or  the  theory  incorporates  a correction  for  refraction. 


3.  EXPERIMEirerAL  FACILITIES:  DEVELOPMEIOT  AND  ARRANGEMENT 


A considerable  portion  of  the  project  has  been  devoted  to  the 
development  of  a Laser  Doppler  Velocimeter  (LDV)  system  that  could  meet  the 
constraints  imposed  by  the  natiure  of  the  experiment: 

(i)  the  system  sho\ald  not  interfere  with  the  flow  field  nor  the 
sound  field, 

(ii)  the  frequency  response  shovild  extend  to  the  vqjper  audio  range, 

(iii)  the  seeding  generator  needed  to  supply  a sufficient  number  of 
scattering  agents  for  the  Laser  Doppler  Velocimeter  must  be 
very  quiet. 

The  development  of  the  latter  proved  to  be  a rather  time-consuming  endeavour. 
3-1  Laser  Doppler  Velocimeter 
3.1.1  Operating  Principle 


A small  particle  travelling  with  a non-relativistic  velocity 
= ^/St  crosses  a light  beam  of  frequency  Vq  and  wave  number  vector  k^, 
1^1  = Ztt  vo/c,  and  scatters  the  spherical  wave  described  by 


-i(2TrVot  - k^-y)  ik|x-yj 


(3.1) 


to  the  observer  at  x in  the  far  field.  S is  a directivity  factor  that  depends 
on  the  particle  siz^,  refractive  index,  and  the  incident  radiation  (Ref.  53); 
k is  the  scattered  wave  number.  The  far- field  approximation  is: 

-i  Jl)  . y 

'sera'*’?’*' ” *t.  ' <3.2) 


where 
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If  a second  beam  of  light  of  the  same  freqviency  Vq  but  with  a different  wave 
vector  is  scattered  simultaneously,  then  the  signal  at  jc  becomes; 
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It  has  been  assumed  that  S varies  sufficiently  slowly  over  the  reinge  0 = 
tan"l(  Iki  - k2|/2|ki  + k2|)  to  be  represented  by  a constant.  The  eaqpression 
Es(jj^  can  be  rewritten  in  the  form  of  an  an^plitude  modulated  signal: 
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The  first  term  is  a wave  with  frequency 
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''  = "o  ^ ^ - 4 


-,p  o’ 


ttrl 


(3.14) 


(3.5) 


and  can  be  thought  of  as  the  high  frequency  carrier  signal.  The  amplitude 
modulation  consists  of  a DC  component  and  a time  varying  signal  with  frequency 


*■0  ■ i?  <4  - 4>  • 5p  <3.6) 

The  latter  is  the  so-called  Doppler  signal;  it  is  actually  the 
difference  of  the  two  Doppler  frequencies  that  are  generated  by  the  scattering 
process.  The  frequency  vd,  which  is  much  lower  than  v,  can  be  detected  by  a 
photomultiplier  and  measured  by  conventional  frequency  analyzers  or  high  speed 
counters.  The  Doppler  frequency; 


'’d 

is  unique  to  the  system  and  varies  linearly  with  the  speed  of  the  particle  in 
the  direction  of  the  ^ ^ vector. 


I'd  -ii  Ife  -41  • iipl  • ♦ 

Vjj  = ^ sin  0 IVp  |cos  0 


(3.7) 
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The  sense  of  direction  is  lost,  hut  can  be  recovered  by  a system  that  uses 
freqviency  shift  techniques. 


For  a more  heuristic  description  consider  the  intersection  of  two 
coherent  light  beams.  At  the  beam  intersection  an  interference  pattern  is 
created  (Fig.  3*1) • The  angle  of  intersection  and  the  wavelength  of  the  light 
determine  the  fringe  spacing: 

d^  = •^  (sin  0)  ^ (3.8) 

A small  particle  crossing  the  fringe  pattern  will  scatter  a modulated  signal, 
the  frequency  of  which  is  deteimlned  by  the  fringe  spacing  and  the  particle 
velocity  normal  to  the  fringes.  The  scattered  light  is  then  collected  on  the 
first  dynode  of  a photomultiplier  tube  (Fig.  3*2).  This  type  of  signal  detec- 
tion is  known  as  the  differential  mode.  A phase-locked  loop  or  a high  speed 
counter  can  be  used  to  estimate  the  Dealer  frequency  ^ich  is  directly  propor- 
tional to  the  particle  speed  normal  to  the  fringe  pattern.  There  is  no 
linearization  required,  as  is  the  case  with  the  hot  wire  anemometer.  Aside 
from  small  scattering  agents,  the  flow  remains  iinaffected  by  the  measiirement . 

The  remote  sensing  eliminates  any  probe-flow  interaction  which  nay  generate 
an  iinwanted,  but  well  correlated,  sound  signal.  A comprehensive  survey  of 
the  operating  principles,  advantages,  and  limitations  of  the  LDV  can  be  found 
in  Refs.  5^  and  55* 

3-1.2  Cptical  Package 

The  optical  arrangement  has  been  built  ip  from  standard  DISA  com- 
ponents consisting  of  a Beam  Splitter,  55l87  Beam  Separator  and  Focusing 

Lens,  and  a 55L12  Photomultiplier,  as  shown  by  the  schematic  representation  in 
Fig.  3*3 • For  minimum  acoustic  interference  a long  range  back-scatter  system 
would  have  been  preferable.  Small  particles  that  follow  the  turbiiLent  flow 
most  faithfully,  scatter  the  incident  beam  in  the  forward  direction;  therefore, 
the  differential  mode  of  light  collection  (Fig.  3 *2)  offers  a better  signal 
strength  (Ref.  56) . 

The  space  available  in  the  anechoic  room  restricted  the  focal  lengths 
of  the  optics  to  be  30  cm.  In  order  to  minimize  acoustical  interference,  all 
optical  components  hswe  been  moved  as  far  away  from  the  flow  as  possible, 
subject  to  spatial  and  vibrational  constraints. 

After  several  design  iterations  a self-contained  optical  package 
was  built.  The  system  can  be  easily  aligned  and  allows  the  beam  intersection 
to  be  placed  virtually  anjnrtiere  in  the  Jet  flow.  There  are  two  main  components: 
the  traversing  gear  and  the  basic  optical  system  (Fig.  3*^)* 

The  X-Y  traversing  gear  provides  a stable  platform  for  the  Laser 
Doppler  optical  system.  The  traversing  gear  has  a range  of  travel  of  30  cm 
along  the  Jet  axis  and  20  cm  laterally  to  it.  An  alimiin\mi  U channel  beam 
svpports  the  beam  splitter  cptics  and  the  photomultiplier  assembly.  The 
horizontal  sipport  is  attached  to  the  traversing  gear  via  a "Lazy  Susan' 
bearing,  providing  the  capability  to  rotate  the  entire  assemibly  about  the 
probe  volume.  This  feature  makes  the  measurement  of  v^  possible.  The 
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sv^Jport  beam  also  serves  as  a platform  for  a Spectra  Physics  12kk  15  mW  He-Ne 
laser  (632.8  nm) . Two  front  surface  mirrors  on  adjustable  mountings  steer  the 
laser  beam  into  the  beam  splitter  optics  75  cm  above  the  laser.  The  beam 
splitter  and  i5)per  steering  mirror  are  attached  to  a bracket  whose  vertical 
position  can  be  adjusted.  A carrier  with  lateral  and  vertical  travel  holds 
the  receiving  optics  and  the  photomultiplier  in  place.  The  support  rests  on 
a 50  cm  optical  bench.  This  arrangement  assures  precise  alignment  of  the 
photomultiplier  with  the  probe  volume.  Four  locating  rings,  epoxied  on  the 
si?)port  beam,  lock  the  optical  bench  in  the  correct  position.  The  optical 
bench  also  serves  as  a counterweight  for  the  beam  splitter  optics,  so  that 
there  is  no  external  torque  on  the  traversixag  gear.  The  entire  system  is 
quite  stable,  and  only  severe  shocks  will  result  in  misalignment.  Prior  to 
installation  in  the  anechoic  room  the  laser  and  the  optical  path  to  the  beam 
splitter  were  covered  with  a cardboard  mantle  to  reduce  dust  accxanulation  on 
the  mirrors.  Figure  3*5  shows  the  optical  system  after  installation  in  the 
anechoic  room.  When  aligned  to  measure  the  velocity  con5)onent  at  4o*  to  the 
Jet  axis,  the  beam  splitter  optics  are  Just  upstream  of  the  nozzle,  thereby 
eliminating  any  interference  with  sound  radiated  in  the  direction  of  the 
microphone . 


3.1.3  Laser  Doppler  Processor 

The  detected  Doppler  signals  are  analyzed  by  a DISA  55L90  counter 
processor.  The  55L90  is  a sophisticated  high  speed  electronic  counter  capable 
of  resolving  5 nsec  intervals  (ICX)  MHz) . The  photomultiplier  signal  is  con- 
ditioned by  adjustable  high  and  low  pass  filters.  The  high  pass  removes  the 
so-called  Doppler  pedestal  which  does  not  contain  any  velocity  information. 

The  low  pass  filter  in^jroves  the  signal  to  noise  ratio  of  the  Doppler  signal 
by  minimizing  the  noise  bandwidth  of  the  photomultiplier  shot  noise. 

A high  speed  Schmitt  trigger  converts  the  filtered  Doppler  signal 
to  a square  wave  that  is  conpatible  with  the  digital  counter  circuitry  used 
to  measure  the  time  ^ taken  by  a single  particle  to  cross  8 fringes  in  the 
probe  volume,  defined  by  the  beam  intersection. 


J The  reciprocal  of  ^ is  8 times  the  Doppler  frequency  which  is 

; directly  proportional  to  the  particle  speed.  The  latter  is  available  at  the 

output  of  a digital  to  analog  converter.  The  DAC  holds  the  old  output  level 
t until  a new  valid  velocity  measxirement  has  been  performed.  A 10  MHz  Doppler 

signal  can  generate  new  output  data  at  the  rate  of  625,000  sauples  per  second. 


In  practice  the  number  of  particles  that  cross  the  probe  volume 
per  unit  time  and  the  orientation  of  the  probe  volume  with  respect  to  the  mean 
flow  direction  govern  the  data  rate  (Appendix  D) . The  main  sources  of  error 
in  the  velocity  estimates  are  large  particles  that  do  not  follow  the  flow,  but 
scatter  a large  amount  of  light,  particles  that  cross  less  than  eight  fringes 
and  multiple  particles  in  the  probe  voliane.  An  internal  validation  system 
minimizes  the  errors  caused  by  these  invalid  signals . For  optimum  data 
acquisition  both  the  data  rate  and  the  percentage  of  valid  data  contained  in 
the  input  signal  must  be  large.  In  this  instance  the  best  use  is  made  of 
particles  that  pass  through  the  probe  voliane.  A digital  subsystem  of  the 
processor  displays  the  output  data  rate  as  well  as  the  validation  percentage. 
Both  these  parameters  were  carefully  monitored  during  each  ea^eriment. 
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3.1.4  Seeding  Generator 

In  order  to  assure  a high  data  rate,  a proportionately  large  nunfcer 
of  scattering  particles  must  cross  the  probe  voluue.  A wide  variety  of 
particles  and  particle  generators  can  be  used  (Ref.  57) . In  this  particular 
application,  the  particle  generator  had  to  be  efficient  as  well  as  quiet.  A 
great  deal  of  effort  was  devoted  to  the  development  of  such  a device,  culminating 
in  two  particle  generators. 

The  first  seeding  generator  uses  cctnmercially  available  .5  alumina 
polishing  powder  as  the  scattering  agents.  The  particles  are  kept  in  a holding 
tank  (Fig.  3.6).  A yanan  air  jet  stirs  the  air  above  the  particles  and  some  of 
them  are  carried  away  with  the  exhaust  flow.  The  point  of  particle  injection 
is  just  ahead  of  the  entrance  to  the  model  jet  plenum.  A fine  stainless  steel 
screen  serves  as  a filter  to  reduce  the  chance  of  particle  coagxalation.  The 
particle  concentration  at  the  nozzle  was  uniform,  and  data  rates  in  excess  of 
60  KHz  were  measured  as  far  downstream  as  7 jet  diameters.  It  was  found, 
however,  that  the  alvonina  particles  would  migrate  into  the  control  room  and 
a^acent  laboratory  space.  The  serious  health  hazard  posed  by  the  fine  alinnina 
dust  prevented  use  of  this  particle  generator. 

By  far  the  safest  seeding  material  is  water,  since  the  water  droplets 
eventually  evaporate.  Yanta  et  al  (Ref.  58),  for  exaii5)le,  used  water  droplets 
to  seed  a s\q)ersonic  flow.  The  most  efficient  method  of  generating  water 
droplets  is  the  process  of  atomization.  Liqiild  is  drawn  into  the  low  pressure 
region  of  a high  speed  jet  flow,  and  the  shear  action  of  the  flow  breaks  off 
smal 3 droplets.  The  size  and  concentration  of  the  particles  is  a function  of 
the  jet  velocity  and  the  position  and  diameter  of  the  liqtdd  simply  line. 

To  assure  a large  particle  concentration,  the  jet  must  possess  a 
large  injector  effect.  This  can  be  achieved  with  a high  speed  jet,  making  it 
necessary  to  find  the  best  ccn^romlse  between  efficiency  and  minimi  noise 
generation.  It  was  found  that  the  nozzle  of  the  atomizing  jet  could  be  very 
small  without  affecting  the  particle  concentration.  The  small  jet  diameter 
assured  that  the  atomizer  would  radiate  a negligible  amount  of  noise;  as  is 
suggested  by  the  U®d2  acoustic  power  law.  The  final  design  is  illustrated  in 
Fig-  3-7.  The  atomizer  consists  of  a miniature  jet  connected  to  a variable 
high  pressure  air  siqjply.  The  liquid  siqjply  line  is  terminated  in  a No.  19 
stainless  steel  hypodermic  tidje  and  positioned  along  the  jet  axis  until  optimum 
seeding  conditions  have  been  reached. 

The  procedure  is  one  of  trial  and  error.  The  air  jet  has  a suffi*- 
ciently  large  injector  effect  to  allow  the  particle  generator  to  be  self  steurting. 
The  seeding  generator  is  placed  in  the  upstream  end  of  the  model  jet  plenum. 

In  static  tests  the  plume  of  the  minijet  spread  to  the  walls  of  the  plenum  in 
a distance  of  30  cm.  This  is  not  expected  to  be  much  different  when  the  model 
jert  is  tvimed  on. 

The  heavier  droplets  settle  more  quickly,  strike  the  walls  of  the 
plenum,  and  are  collected  in  two  drains.  The  atomizer  fluid  is  stored  in  an 
airtight  glass  jar  outside  the  plenum.  A pressure  equalization  line  maximizes 
the  pressure  drop  between  the  reservoir  and  the  miniature  jet.  Based  on  the 
amount  of  water  collected  in  the  drains,  it  is  estimated  that  the  atomizer 
has  an  efficiency  of  70^.  There  is  also  some  loss  due  to  evaporation,  reducing 
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the  efficiency.  Nevertheless,  data  rates  were  fovmd  to  be  consistently  above 
20  KHz  at  all  locations  where  velocity  measioreinents  used  in  the  cross- 
correlations were  performed.  The  overall  water  consuiqption  was  3*3l/hour. 


3.2  Analog  Squaring  Modxile 

3.2.1  Non-Ideal  Behaviour  of  Analog  Multipliers 


An  analog  multiplier  was  used  to  con5)ute  the  signal  proportional  to 
Uk^(X>^) • ^ typical  analog  multiplier  in  the  sqtiaring  mode  is  described  by 

■Uae  equation; 


Sv., 


u 

w 


(3.9) 


where  Vu  is  the  input  voltage  and  Vu2  the  output  voltage.  The  second  term  is 
often  referred  to  as  AC  feed-through.  The  factor  7 is  a function  of  frequency 
and  signal  anplitude. 


Even  high  performance  multipliers  suffer  from  AC  feed-through.  In 
many  instances  the  non- ideal  behaviour  can  be  tolerated  as  a negligible  error. 
This  does  not  apply  to  cross-correlation  experiments.  The  ii^ut  voltage  Vu 
is  proportional  to  the  velocity  %(2^,t) , the  constant  of  proportionality  being 
the  sensitivity  of  the  transducer  system: 


V = k u (y,t) 
u u ' 


(3.10) 


The  cross- correlation  function  of  V^2  pressure  signal  Vp  is: 

V oV  = -7^  k^  u^p  + vkk  -I-  up 
^2  p lOV  upx-^  'updrxr 


(3.11) 


and  the  cross- spectraJ  density  function  is; 


(3.12) 


In  order  to  ass^lre  a good  estimate  of  cpuOp(uj),  the  leading  term  must  dcninate 
at  all  frequencies; 


k 'Pu  2p^“^ 
xf 


(3.13) 


From  theoretical  considerations  it  can  be  shown  that 

\ ZpM 

--  - > -i  = 2U  3P  (u) 


(3.14) 
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where  ^(u)  is  the  ratio  of  the  self  and  shear  noise  spectra.  If  the  inequality 

k U 

<3.15) 


is  not  satisfied,  then  the  estsimate  of  the  cross- spectral  density  cPujj2p  as  well 
as  the  cross-correlation  will  be  in  error. 

3-2.2  Squaring  Circuit 

As  discussed  above,  analog  multipliers  are  not  ideal  elements,  but 
suffer  from  AC  feed-through.  After  testing  several  other  models,  a Burr-Brown 
4206  multiplier  vjas  selected  on  the  basis  of  good  AC  feed-through  rejection. 

In  order  to  assure  that  the  criterion  stated  in  Eq.  3-15  is  always  satisfied, 
the  performance  of  the  multiplier  has  been  iii5)roved.  The  AC  feed-throxigh  can 
be  eliminated  if  a signal  equal  to  -7  S/St  {Vxi)  is  added  to  the  output  of  the 
multiplier.  As  7 is  a function  of  signal  aaqplitude,  one  can  only  hope  for 
limited  success.  The  circuit  diagram  of  the  squaring  module  with  AC  feed- 
throxjgh  correction  is  shown  in  Fig.,3-8-  The  multiplier  has  least  error,  when 
it  is  operated  in  a negative  quadrant  (Vx*Vy  is  always  negative) . The  basic 
squarer  is  built  up  from  a buffer  an^ilifier  Q]  and  a unity  gain  inverter  Q2 
connected  to  the  X and  Y inpaits  of  the  4206  respectively.  A differentiator 
followed  by  a variable  gain  aii5)lifier  (jjo  to  Q5)  contputes  the  error  signal  to 
be  added  to  the  output  of  the  squarer.  Samnning  aii5)lifier  performs  the 
addition.  Switch  Si  allows  the  basic  sq\iarer  to  be  adjusted  independently  of 
the  error  can5)ensatlng  network.  A relative  iinprovement  of  10  dB  is  achieved 
over  a wide  range  of  frequencies,  when  the  system,  is  properly  trimmed  (Fig. 

3.9). 


3.3  Acoustical  System 
3.3.1  Anechoic  Room 

All  measurements  were  performed  in  the  anechoic  room  at  the  University 
of  Toronto  Institute  for  Aerospace  Studies.  The  anechoic  room  has  dimensions 

4.2  X 2.9  X 2.1  m?  between  the  tips  of  the  fibreglass  wedges  which  have  a wedge 
depth  of  20  cm.  A .32  mm  thick  lead  sheet  provides  extra  transmission  loss 
for  external  sound  that  may  be  transmitted  through  the  concrete  walls.  The 
cut-off  frequency  of  the  room  (defined  by  deviation  from  the  l/r  law)  is  approxi- 
mately 300  Hz.  The  background  noise  level  is  less  than  28  dBA. 

A fibreglass  lined  duct  ( .75  x .75  cross-sectional  area)  serves  as 
an  exhaust  for  the  model  jet  installed  in  the  anechoic  room.  All  equipment  is 
wrapped  in  fibreglass,  to  eliminate  excessive  reflections.  A detailed  description 
of  the  anechoic  room  can  be  found  in  Ref.  59- 

3.3.2  Model  Air  Jet 


The  model  air  jet  used  in  the  experiments  is  capable  of  producing 
nozzle  velocities  from  10  m/sec  to  300  m/sec.  The  nozzle  has  an  exit  diameter 
of  1.91  cm  (3/4  in)  and  the  area  ratio  (the  square  of  the  ratio  of  the  plenum 
diameter  most  upstream  of  the  nozzle  and  the  orifice  diameter)  is  110.  The 
nozzle  contovu:  has  been  designed  to  provide  slug  flow  at  the  orifice  (Ref. 

59) . Aim  long  steel  pipe  acts  as  a plentun. 
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The  model  air  Jet  is  connected  to  a fibreglass-lined  muffler  in 
the  control  room  via  a flexible  air  hose.  The  muffler  isolates  the  Jet  from 
an  adjustable  two  stage  pressure  regulator.  A catcpressor  st^lies  70  to 
100  psi  air  to  the  system.  The  ccn^>ressed  air  is  passed  through  a water  trap 
and  regxjlated  down  to  6o  psi  by  the  first  regulator  stage.  The  second  stage 
is  a Fisher  Type  99  Regulator  with  feedback  stabilization.  Cross- correlations 
of  the  valve  noise  of  the  regulators  and  the  Jet  noise  indicate  that  the  muffler 
and  the  flexible  air  line  provide  conplete  acoustic  isolation  from  valve  noise. 

The  static  pressxire  i;pstream  of  the  nozzle  is  monitored  by  a mercxory 
manometer.  All  cross- correlation  e:q)eriments  were  performed  at  a static  head 
of  50  ram  Hg,  corresponding  to  a Jet  velocity  of  105  m/sec. 

3. 3 *3  Microphone  System 

A Bruel  & KJaer  l/2"  condenser  microphone  served  as  the  pressure 
transducer.  The  l/2"  microphone  type  433^  is  the  best  compromise  between 
sensitivity  (12.5  mV/Pa)  and  frequency  response  (4  Hz  to  40  KHz) . The  micro- 
phone is  coupled  to  a Brviel  & KJaer  2619  FEC  preaii5)lifier  driven  by  a Bruel  & 

KJsier  2801  microphone  power  si:i)ply  via  a special  low  capacitance  10  m extension 
cable.  A 2.5  m long  balanced  aluminum  T beam  supports  the  microphone  and 
preanplifier  assembly.  The  microphone  boom  can  be  rotated  about  its  vertical 
axis  which  is  aligned  with  the  Jet  centre  line  and  the  nozzle  exit  plane.  The 
microphone  is  located  in  the  horizontal  plane  of  the  Jet  at  a distance  of  I.60  m 
from  the  Jet  nozzle  and  at  an  angle  of  40*  to  the  Jet  axis,  when  cross-correlation 
experiments  ane  performed.  The  average  time  delay  between  the  time  of  emission 
of  the  signal  in  the  Jet  flow  and  reception  at  the  microphone  is  4.5  msec. 

Careful  attention  has  been  paid  to  eliminate  ground  loops  caused  by 
possible  electrical  contact  of  the  microphone  system  with  its  metal  svpports. 

A Keithley  Model  102  B decade  isolation  anplifier  boosts  the  low  level  micro- 
phone signal  by  a factor  of  80.  The  l/2"  microphone  limits  the  overall 
frequency  response  to  be  4 Hz  to  4o  KHz. 

3.4  Hot  Wire  Anemometer 


A Thermosystems  TSI  1034  A lineeirized  anemometer  was  used  to  measure 
the  tvorbulence  and  the  mean  flow  characteristics  of  the  model  air  Jet.  These 
measiurements  served  as  reference  data.  The  anemcmeter  has  a frequency  response 
from  DC  to  16  KHz,  measured  by  the  sine  wave  method  (Ref.  60) . All  measurements 
were  performed  with  hot  film  probes  type  1210  AG  - 10  featuring  a .0254  mm 
sensor  diameter  and  a 1.27  inm  sensor  length.  These  particular  probes  have  been 
developed  by  Lee  (Ref.  6I)  and  reduce  the  probe  flow  interference  found  in 
conventional  short  stem  hot  film  probes. 

The  water  trap  of  the  air  supply  to  the  model  Jet  does  not  remove 
very  .■gtiai i oil  droplets  that  are  picked  ip  by  the  air  in  the  coapressor  stage. 
The  particle  concentration  is  not  significant;  however,  when  a hot  film  probe 
is  exposed  to  the  Jet  flow  for  an  extended  period  of  tine,  it  is  covered  by 
a baked-on  oil  film.  After  use,  the  sensors  were  cleaned  with  a fine  brush 
dipped  in  Veursol;  the  solvent  cocibined  with  a gentle  brushing  action  removed 
the  oil  film.  The  hot  film  probes  are  surprisingly  riigged  and  can  stand  ip 
to  r^eated  vigorous  cleansing. 
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3*5  Signal  Processing  Instrumentation 


Two  principal  coaoi^utatlon  instruments  were  used  dvuring  the  inves- 
tigation; a Princeton  Applied  Research  101  Correlation  Function  Cco5)uter  and 
a Spectral  Dynamics  DSP  36O  sigpial  processor. 

3.5.1  PAR  101  Correlation  Function  Coiqputer 

The  PAR  101  is  a hybrid  (part  analog,  part  digital)  instrument  that 
cc(n5)utes  a naming  time  average  of  100  lag  products.  The  result  is  an  estimate 
of  the  correlation  function  at  100  eqiially  spaced  points  in  the  time  span 
T = To  to  T = To  + T.  The  precaii5)utation  delay  To  can  be  varied  in  integral 
steps  of  T,  the  delay  range.  The  integration  time  of  the  RC  averagers  is 
adjustable.  In  order  to  be  able  to  measure  correlations  ais  low  as  .005,  it 
is  necessary  to  \xse  a 400  sec  integration  time  and  wait  for  at  least  4 time 
constants  before  saii?)ling  the  output.  The  inherent  noise  floor  of  the  instru- 
ment makes  it  iii?)ossible  to  detect  correlations  with  correlation  coefficients 
less  than  .002. 

3-5.2  Spectral  Dynamics  Digital  Signal  Processor  36O  (DSP  360) 

The  DSP  360  is  a two  channel  fast  Fourier  transform  computer  capable 
of  performing  cross-spectral  density  and  cross- correlation  fiaiction  estimates. 
This  particular  instrument  became  available  in  the  later  stages  of  the  investi- 
gation. All  Laser  Doppler  Velocimeter  based  measurements  were  analyzed  with 
the  DSP  360. 

The  fast  Fourier  transform  cooa^juter  estimates  the  cross- spectral 
density  directly  (Appendix  C) , resulting  in  a significant  decrease  of  con5)u- 
tation  time,  when  coii5)ared  to  the  method  eii5)loyed  by  Lee  (Ref.  61) . Memory 
limitations  restrict  the  nimber  of  averages  that  can  be  performed  on  the 
output  data  to  4096,  equivalent  to  a 100  sec  averaging  time,  when  a 0 to  12 
KHz  analysis  range  is  selected. 

Despite  the  averaging,  the  cross-spectral  density  functions  were 
not  free  from  noise.  It  was  formd  that  the  cross- spectral  density  estimates 
could  be  aneilyzed  with  least  error,  when  they  were  displayed  in  an  anmlitude- 
phase  angle  format.  Such  a display  allows  direct  multiplication  by  f2,  if 
the  function  is  plotted  against  logarithmic  frequency. 

For  computational  reasons,  the  analyzer  uses  only  half  the  amount 
of  data,  \rtien  estimates  of  correlations  are  being  0005) uted.  The  decrease  is 
effective  averaging  time  brings  about  a poorer  signal  to  noise  ratio  for  the 
cross -correlation  functions. 

3.5*3  Ancillary  Instrianentation 

In  order  to  be  able  to  normalize  the  correlation  fvinctlons,  the 
RIG  voltages  of  the  input  signals  must  be  known.  The  input  voltages  to  the 
correlation  computers  were  measured  with  two  Bruel  & Kjaer  2417  Random  Noise 
Voltmeters,  each  with  a 100  sec  integration  time.  The  RIG  voltages  eilso 
served  as  a check  on  possible  drift  in  the  measiirlng  Instrumentation.  The 
coo5)Uted  correlations  and  spectra  were  drawn  by  a Mosely  2D4  X-Y  recorder. 

The  hard-copy  output  data  formed  the  basis  for  further  data  anEiLysis. 
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4.  BgPERIMEagAL  CONSIDEEtATIONS  AMD  CROSS-CHECKS 


Before  attenqpting  to  measiire  the  cross-correlations  and  cross- 
spectral  densities,  it  is  necessary  to  accovint  for  mechanisms,  other  than 
described  by  the  theory,  that  are  capable  of  affecting  the  cross-correlations/ 
spectra.  Some  of  these  can  be  minimized,  others,  hopefully,  are  not  signifi- 
cant. The  main  sources  of  error  that  are  likely  to  occur  during  the 
estimation  of  the  functions  Ru^(x,y,T),  Ru3f2p(x,y,T) , cp^^^p(x,y,u)) , and 
9ux2p(x,^,u))  are;  xP--k 

(i)  distortion  of  signals  by  the  transducers  and  an^Jlifiers, 

(ii)  distortion  of  sound  waves  by  propagation  through  a turbulent 
shear  flow, 

(iil)  spiLTious  sound  generated  by  the  seeded  flow. 

4.1  The  Effect  of  Transducers  and  ABS)liflers 

The  can5)utatlon  of  the  cross- correlation  functions  and  the  cross- 
spectral  densities  is  performed  with  electrical  signals  (voltages)  that  are 
near  replicates  of  the  physical  quantities  they  represent.  Ideally  a trans- 
ducer should  convert  a physical  property  of  a medium  into  an  electrical  signal 
for  all  possible  states  of  the  particular  property  being  measured.  This 
requirement  dictates  that  the  transducer  have  a flat  response  over  all  frequen- 
cies, an  infinite  dynamic  range,  and  be  free  of  noise. 

A practical  coapromise  is  to  use  transducers  that  have  a smooth 
frequency  response  at  least  one  decade  beyond  the  lowest  and  highest  frequencies 
of  interest  to  the  experimenter.  The  transducers  must  be  operated  well  within 
their  dynamic  range,  in  order  to  minimize  electrical  noise  and  signal  clipping. 
Similar  considerations  apply  to  anplifiers. 

The  action  of  a transducer  is  analogous  to  passing  the  signal  to  be 
measured  through  an  electrical  network.  For  a linear  transducer,  the  electrical 
network  is  made  up  of  passive  conponents  such  as  resistors,  capacitors,  and 
inductors;  all  anplifiers  are  assumed  to  be  ideal.  The  output  of  such  a 
transducer  is  given  by; 


g(t) 


t - a)  da 


(4.1) 


•vdiere  f(t)  is  the  physical  quantity  being  measured,  and  h(t)  is  the  inpilse 
response  of  the  equivalent  electrical  network  of  the  transducer  (Ref.  62).  For 
an  ideal  transducer  h(t)  = A 6(t) . 


The  cross -correlation  of  the  output  signeJ.s  of  two  dissimilar  trans- 
ducers measuring  fi(t)  and  f2(t)  respectively  is; 


g^lt)  g2(t  - t)  = J h^(a)  f^(t  - a)  h*(p)  f|(t  - T - p)  da  d^ 

= J h^(a)  hj(p)  f^(”t)  f^(t  - T - p + d)  <»  dP  (4.2) 
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The  correlation  function  of  the  process  ?i(t)  f2(t  - t)  is  distorted  by  the 
action  of  a double  convolution.  The  degree  of  distortion  depends  cm  the 
nature  of  hi(t)  and  h2(t) . 

The  cross- spectral  density  on  the  other  hand  takes  the  form: 


‘PglggC")  = J giU;  ggCt  - t)  e^""  dr 

= J h^(a)hJO)  f;;(t)  fg(t-T-p^-a)  e^"(“-P)dTdadP 

= <4.3) 


The  cross- spectral  density  (Pf^foCw)  process  is  modified  by  the  product 

of  two  transfer  functions  are  the  Fourier  transforms  of  the  respective 

transducer  responses.  A good  transducer  will  be  'flat'  over  some  bandwidth 
of  B Hz,  and  the  anqplitude  of  the  measured  cross- spectrum,  i.e. 

will  be  a good  estimate  of  | over  the  frequency  range  where  | 

is  nearly  constant. 

The  effect  of  severely  mismatched  and  slightly  mismatched  transducers 
is  illustrated  in  Figs.  U.l  and  4.2.  The  case  of  severe  transducer  mismatch  is 
modelled  by  a 'C  weighting  network  and  a I/3  octave  filter  centred  at  1.6  KHz. 
The  signals  to  be  measured  are  generated  by  a cannon  random  noise  source.  The 
cross-correlation  of  the  process  has  two  interesting  properties;  the  correla- 
tion function  'rings'  at  the  fiJter  centre  frequency,  and  the  envelope  of  the 
measured  correlation  peaks  near  .3  msec  time  delay;  the  cross- correlation  that 
would  be  measiured  if  the  transducers  were  ideal  is  symmetric  and  has  a Tn»^v^^n1lnl 
at  zero  time  delay . Even  a slightly  mismatched  transducer  pair,  6is  represented 

by  'C  and  'B'  filters  will  cause  some  distortion  of  the  correlaticm  function 
(Fig.  4.2) . 

The  cross- spectral  density  estimates  can  be  Interpreted  much  more 
readily,  since  the  frequency  response  of  the  transducers  sla5)ly  multiply  the 
cross-spectral  density  to  be  measured.  It  is  more  difficult  to  account  for 
the  effect  of  the  double  convolution.  For  this  reason  it  is  preferable  to 
analyze  data  in  the  frequency  domain.  An  added  benefit  is  that  n fold  differen- 
tiation with  respect  to  time  is  equivalent  to  multiplication  by  (-iu)n  in  the 
frequency  domain. 

When  the  signals  to  be  ancQ.yzed  are  broadband  in  nature,  the  range  of 
time  delay  for  which  the  correlation  is  significantly  greater  than  zero  is  of 
the  order  of  l/B,  where  B is  the  signal  bandwidth.  For  mnim  correlation 
coefficients  the  iii5)ortant  details  of  the  correlation  functions  will  be  buried 
in  instrument  noise,  resulting  in  a poor  estimate.  The  cross- spectral  density 
function  of  the  sa  e process  is  not  immune  from  noise,  but  the  estimate  of  the 
cross- spectral  densj-ty  is  much  less  sensitive  to  noise,  since  it  extends  over 
a wide  range  of  B Hz. 
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4 .2  The  Effect  of  Soiand-Flow  Interaction 


The  sound  generated  by  the  tiirbxilent  fldw  must  travel  through 
turbvilence  as  well  as  a mean  velocity  greidient  before  reaching  still  air. 

The  mean  velocity  gradient  refracts  the  sound  away  from  the  jet  axis  (Ref. 

13) . If  the  field  position  of  the  observer  is  more  than  4o*  from  the  jet 
axis,  then  the  effect  of  refraction  is  minimal.  Scattering  of  sound  by 

tiurb\alent  flow  is  always  present.  As  a resxolt  of  the  interaction  with 
turbulence,  the  sound  wave  emitted  by  a source  in  the  jet  flow  will  suffer 
a Doppler  shift.  The  shift  is  random  because  of  the  randomness  of  the  tur- 
bialence  in  space  and  time.  The  overall  effect  will  be  a spectral  broadening. 
The  analysis  of  the  scattering  problem  is  rather  caii5)lex  (Refs.  63,  64,  65), 
and  it  is  more  expedient  to  perform  a controlled  experiment  that  allows  one 
to  assess  whether  or  not  refraction  and  tvirbiolent  scattering  are  significant 
so\irces  of  error  in  the  estimation  of  the  cross-correlation  and  cross- spectral 
density  fimctions. 

The  experimental  arrangement  of  the  test  is  similar  to  the  one 
employed  by  Atvars  et  al  (Ref.  I3)  in  the  study  of  refraction  of  sound  by  a 
jet  flow  (Fig.  4.3).  An  omnidirectional  point  soiirce  was  placed  in  the  jet 
flow  field  and  driven  by  a pure  tone  signal.  The  sound  received  by  a micro- 
phone in  the  far-field  was  then  conpared  to  the  input  signal.  The  relative 
amplitude  and  phase  of  the  input  and  received  signals  were  measured  for 
several  jet  velocities  (67,  94,  132  m/sec),  source  positions  (yi/D  = 2,  6; 
r/R  = 1)  and  microphone  positions  (-20®  to  90*).  The  microphone  signal  was 
filtered  by  two  narrow  band  frequency  analyzers  that  were  tri  mmed  to  zero  phase 
error  at  the  signal  frequency.  A PAR  129  lock-in  amplifier  served  as  a phase 
detector  and  R^  meter.  The  bandwidth  of  the  detection  system  was  measured  to 
be  .1  Hz. 


The  RMS  output  voltage  of  the  lock-in  amplifier  is  a measure  of  the 
peak  amplitude  of  the  correlation  between  the  signal  emitted  by  the  point 
source  and  the  sound  received  by  the  observer.  The  narrow  bandwidth  of  the 
detection  system  effectively  eliminates  any  uncorrelated  noise  signals. 

Changes  in  the  RMS  output  can  be  caused  by  a combination  of  an  altered 
directivity  pattern  (refraction)  and  spectral  line  broadening  (turbialent 
scattering) . It  can  be  argued  that  the  amount  of  turbulent  scattering  is 
governed  by  the  ratio  of  the  mean  turbulence  scale  and  the  wavelength  of  the 
incident  sound.  The  ratio  is  proportional  to  D/X;  D is  the  jet  diameter. 
Values  of  D/X  at  the  three  test  frequencies  (.069  at  1250  Hz;  .139  at  2500  Hz; 
.278  at  5000  Hz)  suggest  that  the  high  frequency  soxind  will  be  subject  to  more 
intense  t\irbiiLent  scattering.  Refraction  will  resxilt  in  a redistribution  of 
acoustic  energy  in  space,  but  the  frequencies  will  be  preserved. 


The  propagation  of  sound  through  the  tiirbvilent  shear  flow  may  also 
be  accompanied  by  a change  in  the  travel  time  of  the  signal  which  is  equivalent 
to  a change  in  the  relative  phase  of  the  signals.  The  phase  measurements  are 
very  sensitive  to  small  temperature  fluctuations  T = Tanib  + ^ aiid  frequency 
drift  f - fo  + Af: 
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The  reference  test  case  (no  flow)  was  repeated  after  each  parameter  change 
tp  minimize  the  error  in  the  phase  measurements.  Frequency  drift  was  held 
to  .1  Hz/min  by  the  use  of  a stable  oscillator. 

The  nieasurements  are  relatively  insensitive  to  changes  in  the 
position  of  the  source  (this  was  also  found  by  Atvars  et  al)  . The  change  in 
arrival  time  is  governed  by  the  microphone  position,  jet  speed,  and  source 
freqxiency  (Fig.  4.4).  At  4o®  to  the  jet  axis,  the  change  in  arrival  time  is 
less  than  ,02  msec.  Such  a shift  in  time  delay  cannot  be  detected  if  the  basic 
delay  increment  of  the  correlator  is  .08  msec,  as  is  the  case  when  the  cross- 
correlation functions  Ruxp(X}^>''')  and  Rux2p(x,j^,T)  are  being  measured. 

AB5)litude  changes  are  also  functions  of  microphone  position,  jet 
velocity,  and  source  frequency  (Fig.  4.5).  As  expected  from  the  early 
refraction  e:q)eriments  and  the  scaling  outlined  above,  the  low  frequency 
conponents  (1250  Hz)  exhibit  little,  if  any,  measurable  change  in  aaplitude. 
Mid-frequencies  (2500  Hz)  appear  to  be  refracted  at  gmali  angles  to  the  jet 
axis,  with  sign^  recovery  at  larger  angles.  The  increase  in  the  Rie  level 
is  small,  as  the  refracted  energy  is  spread  over  a much  greater  area.  High 
frequencies  (5000  Hz)  suffer  the  largest  aii?)litude  changes.  There  is  no 
ccnplete  recovery,  even  at  90®  to  the  jet  axis.  Apparently  the  turbulent 
scattering  has  broadened  the  spectrum  of  the  pure  tone  to  such  an  extent 
that  it  contains  significant  energy  outside  the  effective  narrow  band  filter 
made  up  of  the  two  ganged  filters  and  the  lock-in  amplifier. 

Even  though  the  line  broadening  is  significant  for  the  high 
frequencies,  it  cannot  be  detected  by  the  fast  Fourier  transform  analyzer 
which  is  operated  at  a frequency  spacing  of  24  Hz,  when  cross- spectral  densities 
or  cross-correlations  are  computed.  The  effective  frequency  spacii\g  of  the  PAR 
correlator  is  of  similar  magnitude.  The  results  suggest  that  the  sound  propa- 
gation through  the  turbulent  jet  flow  will  not  affect  the  estimates  of  the 
cross-correlations  and  cross-spectra  for  the  frequency  range  of  interest. 

4.3  Jet  Flow  Field 

A series  of  measurements  was  performed  with  the  hot  wire  anemometer 
to  determine  the  character  of  the  jet  flow,  and  to  identify  any  possible 
discrepancies  from  typical  jet  flow  behaviour.  The  parameters  of  interest 
are  the  mean  flow  distribution,  the  ttorbulent  flow  distribution,  and  the 
turbulence  spectrum.  Figure  4.6  illustrates  the  mean  velocity  profiles 
Ui(yi,r)/Uj  and  the  turbulence  intensity  profiles  Vui^(y2_,x)/Uj  as  a function 
of  distance  yi  from  the  jet  nozzle.  The  potential  core  of  the*^  jet  extends  to 
four  diameters  downstream  of  the  nozzle,  and  the  mean  velocity  profiles  in 
the  transition  region  e xhibit  the  typical  Gaussian  shape  (Ref.  66). 

In  the  mixing  and  transition  regions  the  radial  profiles  of  turbiilence 
intensity  peak  near  the  r/R  = 1 position.  The  high  turbixlence  levels  gradually 
diffuse  from  a narrow  ridge  centred  about  r/R  = 1 and  'fill  in'  the  region  of 
low  tiirbulence  near  the  jet  axis,  as  is  shown  by  contours  of  equal  turbulence 
intensity  (Fig.  4.7). 

A typical  turbulence  spectrum  (Fig.  4.8)  is  flat  to  some  criticsd 
frequency  f(yi,r)  and  then  decays  at  a rate  of  about  6 dB/octave.  No  measiare- 
ments  were  performed  for  axial  positions  greater  than  yi/D=7.  The  jet  flow 
was  fo;md  to  exhibit  the  characteristic  features  of  a 'typical'  jet  (Ref.  67). 
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4.4  Jet  Noise 


Nossier  and  Ribner  (Ref.  37)  have  devised  a method  which  allows 
one  to  extract  the  self  and  shear  noise  spectra,  as  defined  in  the  Ribner 
theory,  directly  from  far  field  jet  noise  spectra.  The  procedure  relies  on 
the  different  directivity  patterns  and  spectral  contributions  of  the  self  and 
shear  noise . 

Consider  the  siitgplified  self  and  shear  noise  model  with  mean  squeire 
shear  noise  sound  pressure  proportional  to; 


The  two  point  velocity  correlation  differs  from  zero  only  in  a correlation 
volinae  l3.  The  shear  noise  will  scale  as: 


m “x'  “x'  4 
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where  fsjj  is  a characteristic  shear  noise  frequency.  Ux^  may  be  approximated 
as  cos2e. 


Similarly  the  mean  square  self  noise 


(yg-yi)  -x  c'^  |x 


(4 -.7) 


scales  as 
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A spectral  component  of  the  shear  noise  will  vqjon  squaring 

become  e^^^.  Hence  as  a first  approximation:  fsp  = 2fgjj;  for  a random 
process  the  frequency  shift  depends  on  the  statistics  of  the  process  and  may 
not  be  equal  to  2. 

The  sheax  noise,  therefore,  is  expected  to  exhibit  a dipole-lihe 
(cos^e)  directivity  and  dominant  frequency  fsH>  8^*^  ^he  self  noise  shovild  be 
omnidirectional  with  dominant  freqviency  2f(gjj.  At  90*  to  the  jet  axis,  there 
is  no  she£ir  noise  contribution,  allowing  one  to  estimate  the  self  noise. 

A more  rigorous  analysis  due  to  Ribner  (Ref.  10)  is  based  on  a two 
point  correlation  model  with  postiilated  statistics  and  accounts  for  source 
motion.  The  central  results  of  his  walysis  are: 
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(4.9) 


} =■“ 

C = convection  factor  [(l  - M cos  )^  + a^M 

S = Strouhal  nunber  fD/Uj^j, 

(p  = power  spectral  density 

The  overall  Jet  noise  spectiana  measured  at  an  angle  0 to  the  jet  axis  (outside 
the  zone  of  refraction)  is  the  sum  of  a self  pid  shear  noise  spectrm.  The 
effect  of  convection  is  a downstream  bias  (C~4)  and  a Doppler  shift  (CS). 
Furthermore,  the  model  predicts  that  the  self  and  shear  noise  ^ectra  are 
similair  but  with  the  self  noise  shifted  by  one  octave 


The  an^jlitude  factor  p is  predicted  to  be  2. 

If  the  self  and  shear  noise  theory  is  not  a valid  description  of  the 
jet  noise,  then  any  attenpt  to  cast  the  jet  noise  spectra  into  the  form  of 
Eq.  4*9  woiild  resvilt  in  a non-unique  shear  noise  spectrm.  Nossier  and  Ribner 
(Ref.  37)  have  tested  the  self  and  shear  noise  model  against  experimental  jet 

noise  data  and  have  found  that  the  shear  noise  spectriun  is  xmique  and  conforms 

to  the  notion  of  Eq.  4.10,  namely  that  the  self  and  sheax  noise  spectra  match 
and  exhibit  nearly  an  octave  shift.  These  results  sipport  the  major  predictions 
of  the  theory  and  lend  credibility  to  the  existence  of  shear  noise. 

The  basic  self  noise  spectrtnn  can  be  estimated  from  jet  noise  measure- 
ments at  90°  to  the  jet  a^s.  At  angles  other  than  90°,  the  self  noise  is 

changed  by  convection  [C"4(p(cs)]  in  a predictable  manner.  Thus  the  shear  noise 
spectrum  can  be  extracted  by  an  algorithm  formulated  by  Nossier  and  Ribner  (Ref. 
37): 


(CS) 

P9O0 


‘PSHEAR^^®^ 


2 2 
cos  0(1  + cos  0) 


<Pp  (CS)  - 0 . 
0 


(4.11) 


Figvire  4.9  shows  plots  of  cPpQ®  for  three  different  far- field  micro- 
phone positions  at  a constant  jet  velocity  of  10^  m/sec.  The  spectral  peak 
moves  toward  the  lower  frequencies  for  decreasing  angles  with  the  jet  axis 
(reverse  Doppler  shift).  The  measurement  at  0 “90°  contains  some  error,  as 
the  mlcrpphone  was  Just  1^  cm  away  from  the  nearest  fibreglass  wedge,  bringing 
about  some  low  frequency  enhancement.  The  over- estimation  of  the  low  frequency 
cooponents  of  the  self  noise  leads  to  a reduced  shear  noise  estimate. 
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The  self  and  shear  noise  spectra,  calculated  on  the  basis  of  Eq. 

4.11,  eire  plotted  in  Fig.  4.10.  The  ratio  of  the  peak  spectral  aaqplitudes 
of  the  self  and  shear  noise  is  1.71,  and  the  ratio  of  the  peak  frequencies 
is  1.36.  Nossier  and  Ribner  (Ref.  37)  have  found  that  a large  body  of  Jet 
noise  data  exhibits  a one  octave  shift  between  the  self  and  shear  noise  spectral 
peaks.  A 7%  Increase  in  the  peak  an^litude  of  the  shear  noise,  for  example, 
would  shift  the  spectral  peaks  sufficiently  far  apart,  and  the  peak  frequency 
ratio  would  be  closer  to  2.  Even  without  these  hypothetical  corrections,  the 
estimates  ccoapare  with  spectra  computed  by  Nossier  (Ref.  68),  from  piiblished 
Jet  noise  spectra  of  a 1.52"  diameter  Jet  at  Jet  velocity  of  122  m/sec  (Fig. 

4 .11) . The  self  and  shear  noise  spectra  extracted  from  far- field  measurements 
serve  as  reference  data  for  the  self  and  shear  noise  estimates  to  be  cooputed 
from  the  cross- spectral  densities  <Puxp(x>iJ>w)  9u3{2p(x>y>“)  • 

4.5  Laser  Doppler  Velocimeter  Performance  Checks 

The  operating  characteristics  of  the  Laser  Doppler  Velocimeter  (LDV) 
are  described  below.  The  overall  performance  of  the  LDV  is  good,  and  an 
accurate  description  of  the  flxiid  motion  can  be  obtained,  svibject  to  the 
limitations  inherent  in  the  measviring  technique. 

4.5 •!  Beam  Intersection 


After  the  initial  alignment  of  the  optics  was  coopleted,  the  beam 
intersection  was  dotible-checked.  It  is  difficialt  to  detect  small  alignment 
errors,  as  the  beam  diameter  is  about  1 mm.  When  a test  objective  is  placed 
at  the  beam  Intersection,  it  is  possible  to  view  a magnified  image  of  the  probe 
volume  on  an  opaque  screen  (Fig.  4.12).  The  position  of  the  intersecting  beams 
is  adjusted  imtil  the  cross-section  of  the  fringe  pattern  becomes  circular. 
Under  these  conditions  the  test  section  is  an  ellipse  of  revolution  with  major 
axis  10  mm  and  minor  axis  1 mm.  Conpared  to  the  sensing  elements  of  the  hot 
film  probes  used  to  survey  the  Jet  flow  (Sec.  4.3),  the  measuring  volvane 
has  increased  by  a factor  of  30O.  The  \pper  velocity  limit  of  the  LDV  system 
in  the  present  configuration  is  I50  m/sec,  well  above  the  Jet  exit  velocity 
of  105  m/sec. 

4.5.2  Measurement  of  Jet  Velocity 

The  LDV  probe  volume  was  located  at  yi/D  = 2 inside  the  potential 
core  of  the  Jet,  and  the  output  voltage  of  the  processor  was  measxired  as  a 
ftinction  of  the  orientation  of  the  probe  volume  with  respect  to  the  Jet  axis. 
The  normalized  velocity  U(0)/U(O)  is  described  by  the  function  cos9.  There 
is  good  agreement  between  the  measurements  and  the  hypothesis  (Fig.  4.13). 

For  angles  greater  than  60*  data  validation  is  affected  by  the  gmaii  nunfcer 
of  particles  that  cross  at  least  eight  fringes  in  the  probe  volume.  Eventually 
an  orientation  is  reached  (dead  angle)  for  which  no  particle  crosses  eight 
fringes,  preventing  the  processor  from  conputing  a velocity  estimate. 

The  spatial  resolution  of  the  LDV  is  determined  by  the  extent  of  the 
probe  volume.  If  there  is  a significant  change  in  the  mean  flow  velocity  in 
the  spein  of  the  probe  volume,  then  the  DC  output  voltage  will  be  the  average 
mean  velocity  in  the  probe  volume.  [The  determination  of  the  mean  velocity 
in  the  presence  of  a sheeo*  flow  across  the  probe  volume  is  controversial 
(Refs.  69,  70);  however,  if  the  shear  is  small,  then  there  cannot  be  a large 
error  no  matter  which  method  of  mean  flow  estimation  is  used.] 
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A more  seriovis  problem  is  the  creation  of  a spurious  tuibulent 
velocity  cco^ponent.  This  phenomenon,  caused  by  the  san^iling  of  single  particle 
velocities,  is  not  found  in  hot  wire  anemometers,  where  the  heat  transfer  rate 
is  automatically  averaged  by  the  sensor.  In  a turbulent  shear  flow  the  experi- 
menter cannot  identify  whether  the  fluctuating  signal  is  caused  by  an  apparent 
turbulent  velocity  or  a real  one.  The  ambiguity  is  of  no  significance  in  cross- 
correlation  experiments,  as  the  two  velocities  are  real  physical  quantities.  Any 
small  shift  of  the  measiuring  point  within  the  probe  volume  cannot  be  detected. 

The  presence  of  an  apparent  turbulent  velocity  due  to  mean  sheair  will 
increase  the  tiurbulence  intensities  and  broaden  the  region  of  peak  txirbulence 
intensity  that  is  found  near  r/R  = 1 (Fig.  4.l4).  The  LDV  based  mean  velocity 
profiles  agree  with  the  profiles  measured  with  the  hot  wire  anemometer  (Fig. 
4.15).  The  difference  in  the  txirbulence  intensity  profile  can  be  explained  by 
the  presence  of  mean  shear. 

4.5.3  Effect  of  the  LDV  Processor  on  the  Turbulent  Velocity  Spectrum 

The  DISA  55L90  processor  computes  the  speed  of  a single  particle 
and  stores  the  value  in  the  output  memory,  until  a new  valid  speed  has  been 
meastired.  The  output  of  the  processor  resembles  a staircase  with  random  step 
height  and  step  length  (Fig.  4.l6).  The  sampling  theorem  (Ref.  71)  states  that 
no  information  is  lost  if  the  digitization  rate  (here  the  data  rate)  is  at 
least  two  times  the  highest  frequency  of  interest.  The  ip)per  frequency  limit 
was  chosen  to  be  10  KHz,  in  order  to  eillow  for  data  rate  fluctuations  near  the 
boundary  of  the  Jet  flow.  The  peak  frequency  radiated  by  the  jet  is  about  2 KHz, 
at  least  2 octaves  below  the  ipper  frequency  limit.  The  average  data  rate 
observed  during  the  cross-correlation/cross-spectral  density  measurements  was 
greater  than  20  KHz,  the  Nyquist  frequency  for  a 10  KHz  spectral  analysis  range. 
When  the  data  rate  is  artificially  decreased  (by  attenuation  of  the  output  of 
the  photomultiplier) , the  high  frequency  estimate  of  the  turbulence  spectrum 
is  reduced  (Fig.  4.17).  For  the  case  of  high  data  rates  (in  excess  of  20  KHz), 
the  power  spectral  estimates  compare  with  the  ones  computed  from  hot  wire 
anemometer  measurements  (Fig.  4.l8) . The  effect  of  the  data  rate  on  a single 
frequency  component  is  analyzed  in  Appendix  E. 

The  previous  tests  support  the  assertion  of  Cliff  et  al  (Ref.  72), 
namely  that  LDV  and  hot  wire  anemometer  measurements  are  compatible. 

4.5.4  Mass  Addition  Due  to  Seeding  of  the  Flow 

In  order  to  assure  data  rates  that  are  consistently  greater  than 
20  KHz,  water  droplets  must  be  added  to  the  jet  flow  at  a rate  of  .64  gm/sec. 

The  mass  of  the  seeding  particles  increases  the  fluid  density.  As  the  droplets 
have  virtually  no  kinetic  energy  when  they  axe  injected  into  the  plenum,  the 
jet  velocity  must  decrease.  The  droplets  may  also  affect  the  turbulence 
structure  and  can  attenuate  high  frequency  soimd  (Ref.  73)-  Based  on  the 
increase  in  mass,  one  predicts  a 6^  change  in  flow  velocity  and  a 2 dB  change 
in  the  radiated  sound  pressure  level;  the  latter  is  confirmed  by  measiurement . 

The  power  spectral  density  of  the  jet  noise  should  exhibit  a 6^  shift  toward 
the  low  frequencies  (as  a result  of  Strovihal  scaling  fD/Ui),  and  a decrease 
in  spectrvan  level,  vdiich  should  be  uniform,  if  the  absorption  due  to  water 
vapour  and  droplets  is  negligible. 
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There  is  a very  slight  excess  loss  for  frequencies  greater  than 
6.3  KHz  (Fig.  4.19).  Thqre  may  still  be  local  changes  in  the  turbulence 
structTore  of  the  jet  flow.  Any  significant  change  in  the  turbulence  should 
alter  the  near  field  presstu-e  signature  and  the  near-field  far- field  pressure 
cross- correlations.  There  is  no  measurable  change  in  the  normalized  cross- 
correlation coefficients  for  the  seeded  and  unseeded  jet  flow  (Fig.  4.20). 

The  only  effect  cf  any  consequence  is  the  velocity  defect  that  is  caused  by 
the  mass  addition. 


5 . EXPERIMENTAL  PROCEDURE 


5.1  Field  Points 

Even  though  the  contributions  to  the  jet  noise  from  a unit  volume 
of  jet  may  change  with  position  in  the  jet,  the  soxirce  pattern  retains  a top- 
bottom  symnetry.  This  property  allows  one  to  survey  only  field  points  that 
are  located  in  the  \;5>per  two  quadrants'  of  a slice  of  jet,  without  losing  any 
information  about  the  remainder  of. the  slice.  The  field  points  were  located 
on  semicircles  of  radius  R centred  at  yi/D  = n,  n=(3j7)-  The  jet  nozzle 
interfered  with  the  Laser  Doppler  'Velocimeter  optics  for  n < 3-  The  coordinates 
of  the  field  points  were  chosen  to  be  (yi/D,  r/R,  0)  = (n,  1,  mTr/4);  m = (0,  4). 
A more  coir^irehensive  survey  was  performed  at  y]_/D  = 3 and  6,  permitting  a 
better  estimate  of  the  contribution  to  the  far  field  pressure  from  a typical 
slice  of  jet  in  the  mixing  and  transition  regions  of  the  jet  (Fig.  5-1) • 

The  sound  pressure  was  measured  at  a single  position  I60  cm  from 
the  jet  nozzle  and  at  an  angle  of  4o®  to  the  jet  axis.  The  location  offers 
the  best  conroromise  between  the  detectability  of  the  shear  noise  contribution 
(cos^e  + cos4e)  and  errors  due  to  refraction  of  sound  by  the  jet  flow. 

5*2  Signal  Handling 

The  signal  processing  schematic  shown  in  Fig.  5.2  traces  the  paths 
of  the  output  voltages  of  the  LDV  and  microphone  systems.  The  output  of  the 
Laser  Doppler  processor  is  split  into  three  can5)onents.  The  DC  voltage 
Vux  ~ (Ux)  is  measured  by  a calibrated  digital  multimeter.  The  turbulence 
signal,  available  at  the  output  of  an  AC  coupled  amplifier  with  a voltage 
gain  of  8,  or  a voltage  proportional  to  the  square  of  the  instantaneous  tur- 
biilent  velocity,  reaches  the  input  of  the  fast  Fourier  transform  (FFT) 
analyzer  via  a unity  gain  wide  band  filter  (20  Hz  to  40  KHz) . The  microphone 
signal  is  passed  through  a simLleu:  wide  band  filter,  but  with  a voltage  gain 
of  100.  The  cross- correlations  and  the  cross- spectral  densities  are  ocnrputed 
by  the  DSP  36O  FFT  analyzer.  The  input  voltages  to  the  analyzer  are  measured 
by  two  RIC  voltmeters.  The  velocity  and  sovind  signals  are  passed  through 
nearly  identical  networks,  thereby  minimizing  any  spectral  distortion  due  to 
mismatched  frequency  response.  The  voltages  that  represent  the  various 
physical  quantities  to  be  analyzed  are  defined  by: 


Vp  = 8000kp  p(x,t) 

It 

00 

k^  = .040  V/m/sec 

(5-1) 

■ m 

kp  = .0125  V/N/m^ 
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5-3  Correlatictn  Measurements 


The  majority  of  the  cross-correlation  functions  were  con5)Uted  over 
the  delay  range  -12.5  to  12.5  msec.  It  was  possible  to  identify  the  major 
maxima  of  the  cross-correlations,  when  several  realizations  (typically  5) 
were  plotted  on  a common  hard  copy  output.  The  ten?)oral  resolution  of  the 
cross-correlation  fxaictions  could  be  iii?)roved  by  the  xxse  of  a shorter  delay 
range  and  suitable  precomputation  delay.  It  was  fovind,  however,  that  the 
decrease  in  averaging  time  brought  about  unacceptably  large  noise  signals 
that  add  to  the  low  amplitude  cross- correlations.  The  PAR  101  correlator 
which  possesses  an  integration  time  independent  of  the  delay  range  could  not 
improve  the  p<}or  resolirtion. 

5.^  Cross-Spectral  Density  Measturements 

The  cross- spectral  density  estimates  were  affected  by  similar 
signal  to  noise  problems.  The  errors  were  minimized  by  displaying  the 
con^juted  spectra  in  an  amplitude  phase  [log|<p/y)  |;  0®,  <p  = tan"^(lin<p/Re<p)  ] 
format.  The  logarithm  'compressed'  the  noise 'that  added  to  the  signal  and 
allowed  a reasonable  estimate  of  the  amplitude  of  the  cross- spectral  density. 
In  addition,  the  logarithmic  format  enhanced  the  dynamic  range  of  the  output 
memory  of  the  FFT  analyzer  which  is  limited  to  12  bit  words. 

The  analysis  below  shows  that  the  phase  information  is  not  needed 
in  the  estimate  of  the  far  field  power  spectrian,  if  the  two  point  cross- 
correlations  that  are  proportional  to  the  source  strei\gths  are  symmetric. 
Consider  the  auto- correlation  of  the  far-field  shear  noise: 


Rpp  (X’-'-)  ~J  ^ - T + c"^  Ix  - y I)  dy  (5-2) 

SHEAR  V 

The  power  spectrum  is  the  Fourier  transform  of  the  auto- correlation: 
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‘^pp  1^  - y I) 

k-1 ! 
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iu)T  , , 

e dy  dT 


= - <Pu  ® 


dy 


(5.3) 


where  «Pu3{j)(*>y>w)  is  defined  by 


but  u^(y,t)  p(x,t  - t)  has  a maximum  at  t = -c^^  |x  - y|.  The  function 
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f(x,y,T  + c’^lx  - y|)  = U (y,t)  p(x,t  - t) 

-f  -r  -7  x-»  -» 


is  assumed  to  be  symmetric  about  t = -c^  |x  - y|,  and 


(5.5) 


9„  „(y,x,u) 
x^ 


= f(x,y,T  + c"^|x  - y|)  e^"'^  dT  (5-6) 


hj 


= ^(x>y><^)  e 


•iu»c;ilx-j^l 


(5.7) 


(5.8) 


5(x,^,u)  is  wholly  real,  since  f(x,y,T)  is  symmetric  about  T = 0. 

^pp  ~J-  “ {>3^(44>'*')  ® ^ }®  ^ 


SHEAR 


- J'  3(x,y,u)  dy 


(5.9) 


The  FFT  analyzer  measures  an  approximation  to  where, 

theoretically,  (p  = -u)Co“^  - J(|.  The  output  format  thus  allows  in5)lementation 

of  Eq.  5.9.  A similar  analysis  may  be  performed  for  the  case  of  the  self  noise. 

The  measured  cross- spectral  densities  were  processed  by  tracing  a 
smooth  best  fit  curve  on  a translucent  paper,  placed  over  the  hard  copy  drawn 
by  the  X-Y  plotter.  After  applying  the  f^  correction  to  convert  [9uxp(s>^>“) > 
'Pux2p(x,y,u)]  to  ’ 9ux2p(4jyj«)  ]>  cross- spectral  densities 

were  multiplied  by  2VUjj  and  1.25V  respectively  and  plotted  in  a double  logarith- 
mic format.  The  amplitude  scale  is  referenced  to  an  arbitrary  level,  common  to 
all  spectral  measurements. 

The  volume  integral  that  sums  the  contribution  of  all  field  points 
to  the  power  spectral  density  of  the  far- field  pressure  (Eq.  2.19)  can  be 
rearranged  to  be  more  conpatible  with  the  measured  data; 


'Pv  V 
P P 

p 

V = — V[  ^ 

{/  ^ V 1+1-25  |<IV^  V 1 1 f } 

A(yi/D)  “o  ^ p , A J 


(5.10) 


(5.11) 


where  A(y]^/D)  is  the  area  of  a slice  of  Jet  of  thickness  d yi/D  and  A is  the 
area  of  the  Jet  at  y^/D  =3*  If  the  cross- spectral  densities  are  of  the  form 


cp(x,y,f)  = (p(x,yj^,f)  (5.12) 

then 

l<lV^^^^(x,y,f)l  j2Vy^(y)d^(y)  ^ d^ 

^ K pv  I ^iv)  # d ^ (5.13) 

J fo  \ P A " ^ 

where  y is  chosen  to  he  (yi/D,  r/R  =1,  0=0).  The  integrals  over  A(yi/D) 
coEnpute  the  contribution  of  a thin  slice  of  Jet  to  the  total  power  spectral 
densities.  In  particular: 


^ ^(yi/D)  * "" 

“aCr)-/ 

Aiy^/3) 


(5.l4) 


were  estimated  from  detailed  measurements  at  y]_/D  = 3 and.  6.  Values  of 
and  Dp  at  other  Stations  were  obtained  by  interpolation  or  extrapolation.  The 
yi  integration  along  the  Jet  axis  is  replaced  by  a simnnation  with  step  size 
Z!iyi/D  = 1;  thus: 
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(5.15) 


The  above  formalism  is  coonpatible  "with  Eq.  2.21  and  was  used  to  sum  the  contri- 
bution of  the  field  points  in  the  Jet  to  the  power  spectral  density  of  Vp,  which 
is  directly  proportional  to  cpp^p(x,u). 

The  strongest  aliases  in  the  spectra  of  Ux  and  Ux^  are  expected  to  be 
in  the  range  of  10  KHz  and  beyond  (based  on  a 20  KHz  data  rate,  see  Appendix  E); 
therefore,  the  analysis  range  was  chosen  to  span  0 to  10  KHz. 
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6.  RESULTS  MD  DISCUSSION 


6 . 1 Cros  s- Correlations 


The  cross-correlations  Ux(y,t)  pCatjt-rJ  and  p(3c,t-T)  are 

expected  to  peak  near  c^^  •and  decay  rapidly  with  increasing  ? (Appendix 
A) . Measvired  correlations  are  largely  obscured  by  noise,  and  only  the  peak 
values  can  be  determined  with  any  degree  of  accuracy.  The  measurements  are 
repeatable,  eis  is  shown  by  Figs.  6.1  and  6.2.  Each  plot  is  a stqperposition  of 
40  realizations  of  the  'shear'  and  'self  noise  correlations  respectively.  The 
increase  of  the  noise  at  the  extreme  values  of  time  delay  is  an  artifact  of  the 
built-in  amplitude  correction  of  the  DSP  3^0  analyzer. 


The  self  noise  correlations  Rux2p(x>y>T)  *^re  similar  in  their  temporal 
behaviovtr  over  a cross-section  of  Jet,  and  vary  only  in  magnitude  as  the  field 
point  location  is  changed.  The  shear  noise  correlations  Ruxp(^>y>‘’’)  possess  a 
marked  dependence  on  the  position  of  the  source  point  (Fig.  0.3)";  the  behaviour 
illustrated  in  Fig.  6.3  is  typical  for  all  field  points.  Cross-correlations 
measured  with  the  PAR  101  correlation  function  computer,  which  evalirates  the 
cross- correlation  directly,  exhibited  the  same  features. 


The  self  noise  correlations  peak  at  the  expected  time  delay  of 
c3^ [x-^  I , whereas  the  shear  noise  correlations  have  peak  amplitudes  at  about 
.2  msec  earlier  (Fig.  6.4).  The  shear  noise  signal  is  apparently  emitted  by 
the  correlation  volume  prior  to  the  velocity  measurement  at  the  source  point. 
Even  though  the  shear  noise  correlation  is  influenced  by  the  behaviour  of  the 
mean  flow  distribution  (Appendix  A)  , no  model  has  been  fo\jnd  that  can  account 
satisfactorily  for  the  observed  phenomena. 

Shear  noise  correlations  measured  with  a hot  wire  anemometer  (Fig. 
6.5)  do  not  behave  in  a similar  manner  in  that  they  do  peak  at  the  expected 
time  delay  of  Co^|x-^].  This  result  is  compatible  with  the  notion  of  probe- 
flow  interaction  which  appears  to  dominate  the  genuine  flow  noise  generated 
by  the  correlation  voliane  that  surrounds  the  probe  (Appendix  B) . 

The  self  and  shear  noise  correlations  are  strongest  along  the  line 
r/R  = 1,  0 = 0.  The  peak  cross-correlation  coefficients 


Ru^(x,y,T) 

V^(y) 


5^(4) 


(6.1) 


are  the  order  of  .004  (Fig.  6.6) . The  poor  resolution  of  the  cross- correlation 
fvmctions  prevented  any  further  data  processing. 


6.2  Contribution  to  the  Jet  Noise  from  a Unit  Volume  of  Jet 

A set  of  measxired  cross- spectral  densities 

0Yy^^yp(x,^,f)  is  shown  in  Figs.  "6.7  and  6.8.  The  lower  trace  of  each  plot 

is  the  amplitude  of  the  cross- spectral  density  displayed  in  a loglp]  format. 
The  upper  trace  represents  the  phase  angle  in  degrees  (-I80®  to  loO*) . 
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Even  though  both  the  ainplitude  and  the  phase  of  the  measured  spectra 
suffer  from  similar  sources  of  error  (mainly  due  to  finite  averaging  times) , 
one  can  estimate  the  smooth  aii5)litude  of  the  cross- spectral  density  with  some 
confidence.  Repeated  mea3\arements  show  that  the  error  signals  are  random. 

The  linear  display  of  the  phase  angle  is  very  sensitive  to  error  signals  when 
the  phase  angle  of  the  cross- spectral  density  is  near  ± l80*j  and  small  pertur- 
bations cause  a large  amoiant  of  jitter. 


Data  processing,  as  described  in  Section  5-^)  converts  the  cross- 
spectral  densities  to  contributions  to  the  self  and  shear  noise  spectra  from 
a unit  volvime  of  jet  as  shown  in  Figs.  6.9  and  6.10  respectively.  The  spectral 
contributions  and  their  peak  frequencies  are  decreasing  functions  of  the 
distance  from  the  jet  nozzle.  Contrary  to  specxilation , a unit  volume  radiates 
a broad-band  signal;  the  peak  frequency  can  be  considered  as  the  preferred 
frequency.  In  the  range  yq/D  =.3  to  7,  the  peak  frequencies  are  approximated 
by  the  en^iirical  relations: 


W -(5.3-  KHz 

-20^)  KHZ 


(6.2) 


The  average  peak  frequency  ratio  ^ self/^SHEAR  1‘5^.  The  cross- spectral 

density  measurements  indicate  that  the  self  and  shear  noise  contribution  from 
a unit  volume  of  jet  to  the  far  field  jet  noise  at  40°  to  the  jet  axis  are 
approximately  equal. 

6.3  Contribution  to  the  Jet  Noise  from  a Slice  of  Jet 


In  what  follows,  the  terms  source  distributions  and  source  patterns 
are  to  be  interpreted  as  so\irce  power  densities  and  not  as  source  strength 
densities  as  vised  earlier,  e g . on  RHS  of  Eqs  2.4,  2 5,  etc.  The  respective 
self  and  shear  noise  source  strength  densities  are  proportional  to  (52ux/5t2) 
and  Ux(52ux/St2) . On  the  other  hand  contributions  fvo^  a unit  volume  to  the 
fax-field  acoustic  power  forms  (S^Ux^Ux^^T^)  ( S^UxI/xU^Ux/^t^ ) 

which  scale  as  u)^(ux^)^  and  , respectively. 

The  contribution  from  slices  of  jet  at  yq/D  = 3 and  6 have  been 
can5)uted  from  detailed  measurements  at  a grid  of  field  points.  At  a given 
axial  station,  the  spectra  are  similar  in  their  frequency  content  and  differ 
only  in  aa^jlitude  (Fig.  6.II).  After  cross-plotting  the  data,  contours  of 
constant  source  contributions  can  be  obtained.  The  reference  position  in  each 
slice  of  jet  is  located  at  r/R  =1,  0=0°. 

At  yq/D  = 3 the  potential  core  (which  should  contain  no  sovirces) 
is  clearly  defined.  The  source  distributions  of  the  shear  and  self  noise 
(Figs.  6.12,  6.13)  6ire  largely  confined  to  the  quadrants  of  jet  closest  to 
the  observer  and  axe  strongest  along  the  circle  r/R  = 1.  The  source  patterns 
in  the  Proudman  format  are  fixed  to  the  vector  and  move  with  the  observer; 
therefore  the  jet  axisymmetry  will  insure  that  the  mean  square  pressure  is 
axisymmetric.  Because  of  the  top-bottom  symmetry  only  the  upper  portion  of 
the  patterns  axe  illustrated. 
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There  is  a close  affinity  of  the  source  patterns  to  the  functions 
u^l^2  and  (ux^)^.  From  a siiEplistic  argument  these  may  be  considered  to  be 
approximate  indicators  of  the  contributions*  per  unit  volume  to  the  far  field 
sheeir  and  self  noise  power.  Ux^U^^  and  (ux^)2  are  not  axi symmetric,  contrary 
to  the  common  assuD5)tion  in  noise  models.  The  turbulence  is  not  isotropic 
and  the  mean  flow  not  quite  parallel  to  the  jet  axis;  thus  there  will  be  a 
biasing  toward  the  observer  at  x.  The  role  of  the  radial  flow  is  indicated 
in  Fig.  6.15.  The  nominal  source  patterns  do  not  behave  exactly  as  the  two 
point  velocity  correlations  that  govern  the  sound  generation,  and  therefore 
are  not  one  to  one  with  the  measured  source  patterns. 

In  the  transiljion  region  the  properties  of  the  jet  flow  approach 
the  characteristic  behaviour  of  a fully  developed  turbvilent  jet.  The  source 
distributions  are  expected  to  became  more  symmetric  about  the  jet  centre  line. 
At  yi/D  = 6 both  the  self  and  shear  noise  source  patterns  are  strongest  in 
the  quadrants  closest  to  the  observer,  with  a pronounced  bulge  along  the  line 
pointing  in  the  direction  of  the  observer  (Figs.  6.I6,  6.17).  The  mean 
velocity  distribution  Ux  gives  the  shear  noise  pattern  a more  circular 
appearance . As  was  the  case  on  the  mixing  region,  the  nominal  source  patterns 
u^Ux2  and  (u:j^)2  are  similar  to  the  soiirce  distributions  of  the  shear  and  self 
noise  respectively  (Fig.  6.I8). 

6 .4  Jet  Noise  Contribution  from  Field  Points  in  the  Region  3 ■<  Yj^/D  < 7 

The  overall  self  and  shear  noise  spectra  at  40®  to  the  jet  axis  can 
be  constructed  from  contributions  of  radial  slices  of  jet  (see  Section  5-4). 
Figure  6.I9  illiistrates  the  spectral  contributions  from  such  slices  in  the 
range  yi/D  = 3 to  7.  As  the  cross-spectra  were  found  to  differ  only  in  their 
relative  amplitudes  at  a fixed  axial  position  the  curves  are  similar  to  the 
ones  in  Fig.  6.9  and  6.10  but  with  different  magnitudes  to  account  for  the 
overall  emission  from  a slice  of  jet  rather  than  a unit  volinne. 

The  total  self  and  shear  noise  contributions  from  an  individual 
slice,  which  is  directly  proportional  to  the  area  under  the  corresponding 
curves  in  Fig.  6.19,  are  nearly  equal.  The  contributions  to  the  far  field 
jet  noise  are  decreasing  fvinctions  of  axial  position  in  the  range  3<  yi/l>  <.  7 
(Fig.  6.20).  The  hybrid  model  discussed  in  Appendix  A is  shown  for  con5>arison. 
The  model  source  strength  distributions  are  based  on  assumed  frequency  and 
correlation  volume  scaling  and  measured  nominal  source  patterns  €ind  agree 
qualitatively  with  the  normalized  self  and  shear  noise  source  strength  distri- 
butions deduced  from  the  cross- spectral  density  measurements  in  the  region 
sxirveyed  by  the  experiment,  namely  yi/B  = 3 to  7.  This  particular  region 
appears  to  be  responsible  for  about  4o^  of  the  total  jet  noise  heard  by  an 
observer  at  40°  to  the  jet  axis. 


* For  the  model  discussed  in  Appendix  A the  respective  self  and  sheeu:  noise 
contributions  to  the  power  spectr\im  of  the  far  field  jet  noise  are  of  the 
form  (Eqs.  A21  and  P£S)i 


dcPgF  « w^Jl(«jy)  ~ J^(w,y) 

, 2 24-#/  \ 

^SH  ” ^x  ^x  « 


If  jT 1 and  3^2  invariant  across  a slice  of  jet,  then  the  contributions 
to  yie  self/ shear  noise  spectra  axe  governed  by  the  nominal  source  patterns 
(ux^)^  and  Ux^  Ux^  respectively. 


i 
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The  behaviour  of  several  frequency  bands  (bandwidth  2^  Hz)  suggests 
that  the  major  part  of  the  source  distribution  is  confined  to  axial  positions 
less  than  yi/D  = 10  (Fig.  6.21).  Soxirce  location  techniques  such  as  en^iloyed 
by  Grosche  et  al  (Ref.  7^)  and  Fischer  et  al  (Ref.  75)  indicate  that  soxirces 
as  far  downstream  as  15  Jet  diameters  cont.dbute  to  the  Jet  noise.  These 
tests  were,  however,  performed  at  considerably  higher  Mach  nunber  than  the 
present  experiments  (Grosche:  M = 0.7  and  1.0;  Fischer:  M = 0.8).  Grosche's 
results  show  that  the  source  distribution  migrates  towards  the  Jet  nozzle  as 
the  Mach  nimber  is  reduced. 

The  sovLTce  strength  distributions  of  Fisher  et  al  (Ref.  75)  exhibit 
a rather  steep  roll-off  in  the  transition  region  of  the  Jet.  The  rapid  loss 
of  source  strength  is  moderated  by  a slow  asyuptotic  decay  for  yq/D  > 8;  thus 
the  axial  source  strength  distributions  deduced  from  the  measured  cross- spectral 
densities  may  well  be  ccmpatible  with  the  ones  that  would  be  measiored  by  an 
acoustic  soTirce  location  technique.  Acoustic  telescope  measurements  of 
Billingsley  and  Kinns  (Ref.  76)  appear  to  corroborate  the  present  results. 

Their  tests  were  performed  on  a real  (i.e.  hot)  turbojet  exhaust,  and  the 
correspondence  is  thought  to  be  rather  fortuitous . 

The  spectral  contributions  to  the  overall  self  and  shear  noise  are 
estimated  by  direct  summation  over  all  slices  of  jet  that  have  been  surveyed. 

As  was  inferred  above,  only  about  4o^  of  the  overall  Jet  noise  appears  to  be 
accounted  for;  hence,  the  spectra  predicted  from  the  measured  cross- spectral 
densities  are  likely  to  underestimate  the  self  and  shear  noise  spectra  which 
have  been  extracted  from  the  far-field  Jet  noise  via  the  Nossier-Ribner 
technique  (see  Section  4.4).  Based  on  the  behavioior  of  the  sources  contained 
in  the  region  y^/D  = 3 to  7,  it  can  be  argued  that  sources  downstream  of  the 
axial  position  y]_/D  = 7 will  contribute  much  of  their  energy  to  the  low  fre- 
quency parts  of  the  spectra,  and  the  sources  near  the  Jet  nozzle  (yi/D  < 3) 
will  add  to  the  high  frequencies. 

The  self  and  shear  noise  contributions  to  the  overall  Jet  noise 
(in  the  present  investigation  only  the  observer  position  at  40®  to  the  Jet 
axis  is  considered)  can  be  estimated  from  far  field  Jet  noise  data  via  the 
Nossier-Ribner  algorithm  (Ref.  37)  or  from  cross- spectral  densities  as 
measured  herein.  Both  methods  are  consistent  with  Ribner's  self  and  shear  noise 
model  and,  if  the  theory  is  valid,  should  resialt  in  identical  predictions  of  the 
self  and  shear  noise  spectra  subject  to  experimental  constraints.  Should  such 
a closure  be  demonstrated,  then  the  credibility  of  the  theory,  which  is  supported 
by  other  evidence  (for  exaaple:  compatibility  with  the  work  of  Mani,  Ref.  23> 
as  demonstrated  by  Ribner,  Ref.  24,  theoretical  two  point  far  field  pressure 
correlations.  Ref.  28)  is  strengthened. 

The  shear  noise  spectrum  conputed  from  measured  spectral  densities 
and  the  shear  noise  spectrvim  extracted  from  the  far  field  Jet  noise  are 
coupaxed  in  Fig.  6.22.  The  shear  noise  contribution  from  the  limited  volume 
sturveyed  by  the  experiment  is  conparable  to  the  total  radiated  shear  noise. 

There  has  been  no  arbitrary  fit  of  the  spectral  anplitudes;  the  direct 
sunmatlon  described  by  Eq.  5-15  has  been  performed.  It  was  reported  earlier 
that  the  'shear  noise  correlations'  exhibited  a marked  change  over  a radial 
slice  of  Jet;  such  a change  takes  the  form  of  a phase  angle  ^(x,y,to) , as  the 
frequency  distribution  appears  to  be  vinaffected  by  changes  in  the  radial  or 
azimuthal  position  of  the  source  point  (Fig.  6.11) . The  estimate  of  ^ is 
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buried  in  the  phase  angle  computed  by  the  FFT  analyzer.  The  shear  noise 
spectrum 


9, 


PP, 


SHEAR 


e^V'Cx,y,u) 


lx-2,1 

e dy 


(6.3) 


is  approximated  as 


(x,w)  ~ f U (y)  |(p  (x,y,w)  | dy  (6.4) 

P^SHEAR  J u^-r-, 

where  f(x,y,u)  has  been  chosen  to  be  -uCo"^|x-y|.  The  correct  phase  angle  is 
obscured  by  a large  experimental  error;  in  general  V'(x,_jr,u)  ^ 
integrand  in  Eq.  6.4  is  simply  the  absolute  value  of  the  integrand  of  the  more 
exact  shear  noise  spectrum  estimate  given  by  Eq.  6.3.  It  can  be  shown  with 
the  aid  of  the  Schwartz  inequality  that  the  shear  noise  contribution,  as 
computed  by  Eq.  6.4,  overestimates  the  spectrum.  It  is  encouraging,  however, 
that  the  spectrian  estimates  are  of  a comparable  order  of  magnitude.  A longer 
averaging  time  is  likely  to  improve  the  definition  of  the  phase  angle  f(x,y,u) ; 
implementation  of  such  a scheme  reqxiires  a digital  computer,  which  was  not"’ 
readily  available  to  the  investigator. 

The  self  noise  spectrum  computed  from  measured  cross-spectral 
densities  and  the  self  noise  spectrum  extracted  from  the  far  field  j et  noise 
(measuring  the  Jet  noise  spectrum  at  90*  to  the  Jet  axis  and  correcting  for 
the  Doppler  shift  and  convection  amplification  to  obtain  the  self  noise  estimate 
at  40*  to  the  Jet  axis)  are  shown  in  Fig.  6.23.  Tbe  peak  amplitude  of  the 
predicted  self  noise  spectrm  is  found  near  3.15  KHz,  whereas  the  self  noise 
spectrvon  computed  from  the  fsir  field  Jet  noise  via  the  'Nossier-Ribner  technique' 
has  a maximum  neao:  2.2  KHz. 


The  measured  cross- correlations  Ux^(;^,t)  p(2^t-T) , in  contrast  to  the 
shear  noise  correlations  Ruxp(*>i»'’^^ » have  maxima  at  w.e  expected  time  delay  of 
Cq'^I^j-^I;  thus  the  phase  ai^le  introduced  above  can  be  set  to  -wco"^  |x"Xl 
no  significant  error,  and  the  estimate  of  the  self  noise  at  4o*  to  the  Jet  axis 
can  be  treated  with  considerdoly  more  confidence  thiui  the  corresponding  shear 
noise  estimate.  The  energy  content  of  the  predicted  self  noise  spectrum  is 
only  4l^  of  the  overall  self  noise  spectrxan.  This  figure  is  in  good  agreement 
with  the  4^  total  energy  contribution  which  was  inferred  fromi  the  behaviour 
of>»4iie  kxial  source  strength  distribution. 

One  form  of  Ribner's  theory  predicts  that  the  shear  and  self  noise 
spectra  have  the  same  shape  after  an  octave  shift;  the  shift  actually  depends 
on  the  form  of  the  time  correlation  of  the  soiirces  (Refs.  10,  77).  In  Fig. 

6.24  the  shear  noise  spectrvm  extracted  from  the  far  field  noise  has  been 
shifted  in  frequency  and  adjusted  in  amplitude  until  the  spectrum  nearly 
matched  the  self  noise  spectrum.  The  best  fit  was  fo\md  for  a shift  of  1.6. 

The  same  frequency  shift  has  been  applied  to  the  shear  noise  spectrum  predicted 
from  the  cross- spectral  densities.  After  matching  the  peak  amplitudes,  it  is 
found  that  the  two  estimates  of  the  frequency- shifted  shear  noise  spectra  fit 
the  corresponding  self  noise  spectra,  as  predicted  by  the  theory.  Nossier  and 
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Ribner  (Ref.  37)  found  that  a large  body  of  experimental j et  noise  data  on 
average  exhibits  close  to  an  octave  shift . The  frequency  shift  appears  to  be 
a weak  monotonic  function  of  jet  velocity  (Table  I of  Ref.  37),  and  thus  the 
present  result  is  thoiaght  to  be  con5)atible  with  the  one  of  Nossier  and  Ribner 
(Ref.  37). 


For  the  sake  of  coii5)leteness,  the  overall  jet  noise  spectrum  measured 
at  40*  to  the  jet  axis  and  the  sum  of  the  predicted  self  and  shear  noise  con- 
tributions, based  on  the  cross- spectral  density  measurements  performed  herein, 
are  shown  in  Fig.  6.25.  The  apparent  fit  for  frequencies  greater  than  2.5  KHz 
is  for  the  most  part  due  to  the  overestimate  of  the  shear  noise  spectrian. 

6.5  Concluding  Remarks 

Over  the  yeairs  Ribner' s self  and  shear  noise  model,  an  extension  of 
Lighthill's  theory  of  jet  noise,  has  been  demonstrated  to  describe  the  major 
features  of  jet  noise  (Refs.  17,  37,  28,  2h) , but  direct  measurement  of  the 
postvilated  sources  via  a cross- correlation  technique  has  met  only  limited 
success.  The  spxirious  'probe  noise'  generated  by  interaction  of  the  turbulence 
and  the  velocity  sensor  - the  hot  wire  probe  - is  believed  to  be  the  major  source 
of  error.  In  the  present  investigation  the  sources,  which  are  proportional  to 
certeiin  rates  of  momentum  flxix,  have  been  measured  by  a 'remote'  sensing  tech- 
nique made  possible  by  a Laser  Doppler  Velocimeter  that  has  been  adapted  to  the 
special  requirements  of  the  correlation  measurements. 

Cross- correlations  and  cross- spectral  densities  of  the  postulated 
self  and  shear  noise  source  terms  and  the  radiated  sound  have  been  measured. 

For  radiation  to  the  field  point  at  4o*  to  the  jet  axis,  the  effective  instan- 
taneous source  patterns  of  the  self  and  shear  noise,  as  derived  from  the 
cross- spectra,  show  a marked  deviation  from  axisymmetry.  _Nominal  source 
patterns,  proportional  to  Ux^x^  for  the  shear  noise  and  (ux2)2  for  the  self 
noise,  are  qualitatively  similar  to  the  meas\ired  source  patterns.  The  effective 
source  patterns  are  observer  oriented,  and  not  fixed  to  particular  jet  flow 
coordinates;  this  property  assures  that  the  overeQ.1  radiated  mean  square  sound 
pressure  is  axisymmetric,  even  though  the  insteintaneous  sound  pressure  is  not. 

The  self  and  shear  noise  spectra  predicted  from  an  aggregate  of 
cross- spectral  density  measurements  of  the  postulated  source  terms  and  the 
radiated  soxind  are  coii5)atible  with  the  self  and  shear  noise  spectra  extracted 
from  far  field  jet  noise  measTirements.  The  source  region  surveyed  herein 
(yi/D  = 3 to  7)  appears  to  radiate  aboixt  of  the  total  jet  noise  heard  by 
the  observer;  therefore  the  predicted  self  and  shear  noise  spectra  are  not 
one  to  one  with  those  extracted  via  the  'Nossier-Ribner'  algorithm. 

Both  sets  of  self  and  shear  noise  spectra  do  match  very  nearly  in 
shape  and  exhibit  a common  frequency  shift,  and  both  sets  have  co!i5)arable 
absolute  spectral  amplitudes;  this  behaviour  is  predicted  by  the  theory.  The 
con5)atibility  suggests  that  the  theory  is  internally  consistent  and  lends  it 
substantial  further  credibility. 
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FIG.  3.1  MDIRE  PATTEEIN  MDDEL  OF  THE  IHTBRSECTION  OF  TWO  COHEREIW  LIGHT 
BEAJe.  LINES  MAY  BE  THOUGHT  OF  AS  SURFACES  OF  EQUAL  PHASE. 
[REF.  DURST,  F.,  STEVENSON,  W.  H.,  MOIRE  PATTERNS  TO  VISUALLY 
MODEL  LASER  DOPPLER  SIGNALS,  APPLIED  OPEICA,  VOL.  15,  1,  PP. 
137-144  (1976)] 


FIG.  3.2  SCHEMATIC  OF  THE  DIFFERENTIAL  MODE  CF  LASER  DOPPLER  SIGNAL 


DETECTION.  'DOPPLER  BURSTS'  SCATTERED  BY  PARTICLES  TRAVERSING 
THE  BEAM  INTERSECTION  ARE  FOCUSED  ON  THE  PHOTOMULTIPLIER 
DYNCOE. 


BEAM 

SPLITTER 


FLOW  DIRECTION 
BEING  MEASURED 


PHOTOMULTIPLIER 

OPTICS 


15  mW  He-Ne  LASER 


FIG.  3-3  PRnrciPAL  OPTICAL  COMPOIffilJTS  OF  THE  LASER  DOPPLER  VELOCIMETER. 
THE  LASER  BEAM  IS  GUIDED  TO  THE  BEAM  SPLITTER;  DOPPLER 
SIGNALS  SCATTERED  BY  PARTICLES  IN  TliE  BEAM  INTERSECTION 
(PROBE  VOLUME)  ARE  FOCUSED  ON  THE  PHOTOMULTIPLIER. 
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FIG.  3*5  LASER  DOPPLER  VELOCIMEIER  OPTICS  AFTER  INSTALLATION  IN  UTIAS 
ANECHOIC  ROOM. 


.5 jjm  ALUMINA  POWDER 


FIG.  3*6  SEEDING  GEKERATQR  FOR  .5  ALUMINA  POLISHING  POWDER. 


20  PSI  AIR  TO  Hg  MANOMETER 


GENERAL  ARRANGEMENT 


ATOMIZER  DETAIL 


SEEDED 

FLOW 


FIG.  3-7  ATOMIZER  INSTALLED  IN  l/TIAS  3A"  <IEI. 
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FIG.  3-0  ANALOG  SQUARING  CIRCUIT  WITH  AC  FEEDTHROUGH  CORRECTION 


• BASIC  SQUARER 

AC  FEEDTHROUGH  CORRECTION: 
ATYPICAL  PERFORMANCE 

♦ least  error  to  10  KHz 
(SHORT  TERM) 


FIG.  3.9  AC  FEEDTHROUGH  AS  A FUNCTION  OF  FREQUENCY  f . VERTICAL  SCALE: 

RATIO  OF  THE  RJe  OF  THE  RESIDUAL  FRjJQUENCY  AT  f KHz  AND  THE 
SQUARED  SIGNAL  AT  2f  KHz.  NOTE  THE  f+1  DEPENDENCE  OF  THE 
BASIC  SQUARER. 
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FIG.  4.1  M)DEL  EXPERBffiNT  ILLUSTRATIlla  THE  EFFECT  OF  SEVERELY  MISMATCHED 
TRANSDUCERS  ON  THE  ESTIMATES  OF  CROSS- CORRELATION  AND  CROSS- 
SPECTRAL  DENSITY  FUNCTIONS. 


FIG.  4.2  MDDEL  EXPERBCMP  ILLUSTRATING  THE  EFFECT  OF  SLIGHTLY  MISMATCHED 
TRANSDUCERS  ON  THE  ESTIMATES  OF  CROSS -CORRELATION  AND  CROSS- 
SPECTRAL  DENSITY  FUNCTIONS. 
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FIG.  4.3  EXFEIRIMEOTAL  AREIANGEMEKT  FOR  THE  DETECTION  OF  SOUND-FLOW 

INTERACTION.  THE  CHANGE  IN  R^E  LEVEL  AND  PHASE  OF  THE  PURE 
TIME  EMITTED  BY  THE  POINT  SOURCE  IS  MEASURED  FOR  VARIOUS  JET 
VELOCTTIES,  SOURCE  AND  RECEIVER  POSITIONS. 


FIG.  4.4  CHAUGE  IN  THE  TRAVEL  TIME  OF  AN  ACOUSTIC  SIGNAL  EMITTED  BY 
A POINT  SOURCE  IN  A JET  FLOW  TO  A RECEIVER  AS  A FUNCTION  OF 
POSIT  ION  RELATIVE  TO  THE  JET  AXIS . 


FIG.  4.5  NORMALIZED  RJG  LEVEL  OF  AN  ACOUSTIC  SIGNAL  EMITTED  BY  A 

POINT  SOURCE  IN  A JET  FLOW  AS  A FUNCTION  OF  RECEIVER  POSITION 
RELATIVE  TO  THE  JET  AXIS.  REFRACTION  IS  THE  MAJOR  MECHANISM 
THAT  MODIFIES  THE  DIRECTIVITY  OF  THE  RADIATION  EMITTED  BY  THE 
POINT  SOURCE.  (ATVARS  ET  AL,  REF.  13) 


FIG.  4.7 


CONTOURS  OF  CONSTANT  TURBULENCE  INTENSITY 
JET  FLOW. 


IN  THE 


FIG.  4.8  TYPICAL  TURBULENCE  SPECTRUM  OF  t) ; y^/D  = 6,  r/R  = 1. 


FIG.  4.9  FAR  FIELD  JET  NOISE  SPECTRA  OF  THE  UTIAS  3/4"  MODEL  JET  AT 
VARIOUS  OBSERVER  POSITIONS  RELATIVE  TO  THE  JET  AXIS.  JET 
VELOCITY  105  JVSEC,  VERTICAL  SCALE  (LINEAR)  IS  REFERENCED 
TO  AN  ARBITRARY  LEVEL. 


FIG.  4.10  RELATIVE  AMPLITUDES  OF  SELF  AKD  SHEAR  NOISE  SPECTRA 
COMPUTED  FROM  FAR  FIELD  JET  NOISE  MEASUREMENTS 

(FIG.  4.9) . SHEAR  NOISE  CONTRIBUTION  AT  40“ 

TO  THE  JET  AXIS. 


FIG.  4.11  RELATIVE  AMPLITUDE  OF  SELF  AND  SHEAR  NOISE  SPECTRA;  RE- 
PLOTTED FROM  NOSSIER.  [REF.  68,  FIG.  10;  ORIGINAL  FAR  FIELD 
DATA;  AHEYA,  K,  K.,  CORRELATION  AND  PREDICTION  OF  JET  NOISE, 
JOUR.  OF  SOUND  AND  VIBR.,  VOL.  29  (2),  PP.  155-168  (1973); 
JET  DIAMETER  1.52",  JET  VELOCITY  122  KV'SEC] 
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FIG.  4.l4  TURBUUaJCE  miENSITY  PROFILE  OF  THE  JET  FLOW  AT  yi/D  = 6 
CCMPARISON  OF  HOT  WIRE  WITH  LDV. 
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FIG.  4.15  MEAN  VELOCin  PROFILE  CF  THE  JET  FLOW  AT  y^/D  = 6 
CC»!PARISON  OF  HOT  WIRE  WITH  LDV. 
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FIG.  4.l6  TYPICAL  OUTPUT  OF  THE  UN  PROCESSOR.  THE  LAST  VALID 
measurement  IS  DISPLAYED  UNTIL  IT  IS  UPDATED  BY  A NEW 
ONE. 
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FIG.  4.17  EFFECT  OF  THE  LDV  DATA  RATE  ON  THE  ESTIMATE  OF  THE  SPECTRUM 
OF  THE  TURBULENT  FLOW  Ui(y,th  Fl/O  = 6;  r/R  = 0. 


FIG.  4.l8  COMPARISON  OF  SPECTRA  OF  THE  TURBULENT  FLOW  UiCy,t)  MEASURED 
WITH  A HOT  WIRE  ANEMDMETER  AND  AN  LDV.  FIELD  POINT:  yi/D  = 
6,  r/R  = 0. 


FIG.  4.19  FIEID  NOISE  SPECTRA  OF  NORMAL  AND  SEEDED  JET  FLOW. 

the  SEEDED  FLOW  CONTAINS  SMALL  WATER  DROPLETS  (~  nm)  . 
THE  ADDITIONAL  MASS  RESULTS  IN  A VELOCnT  DEFECT  AND 
HENCE  A REDUCTION  IN  SPECTRUM  LEVEL. 
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TIME  DELAY 

FIG.  6.2  MEASURED  CROSS -CQRRELAI IONS  OF  THE  SQUARE  OF  THE  INSTANTANEOUS 
FLOW  v^(y,  t)  AND  THE  FAR  FIELD  JET  NOISE  p(x,  t)  . SOURCE 
POSITION  aa/D  = 6,  r/R  = 0. 
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FIG.  6.3  ILLUSTRATION  OF  THE  BEHAVIOUR  OF  THE  'SHEAR  NOISE  CORRELATION' 
Ruxp(x>  t)  at  yi/D  = 6. 


FK.  6.4  TIME  DELAY  OF  THE  PEAK  VALUES  OF  THE  CROSS-CORRELATION 

T ) (SHEIAR)  and  I^v^2p(*>  T ) (SEU)  . 


FIG.  6.5  NORMALIZED  'SHEAR  NOISE  CORRELATION’  Ruxp(x,  y,  T)//^^(y)pW 
MEASURED  WITH  A HOT  WIRE  ANEMOMETER.  SOUR'3e  POSITION;  yV/D  = tr, 
r/R  = 1,  0 = 0. 


FIG.  6.6  NORMALIZED  CROSS- CORRELATION  COEFFICIEHTS  AS  A FUNCTION  OF 
AXIAL  POSITION.  ALL  SOURCE  POINTS  ARE  LOCATED  ON  THE  LOCUS 
r/R  = 1,  0 = 0“. 
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?IG.  6.7  MEASURED  CROSS-SPECTRAL  DENSITY  OF  THE  TURBULENT  FLOW 
uac(y,  t)  AND  p(x,  t) . SOURCE  POIUr  yi/D  = 6,  r/R  = 1, 
<t>  = 135*;  OBSERm  IS  AT  4o«  TO  THE  JET  AXIS. 
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FIG.  6.8  MEASURiB  CROSS- SPECTRAL  DENSITY  OF  THE  SQUARE  OF  THE  TURBULENT 
PLOW  MxHy,  t)  AND  p(x,  t)  . SOURCE  POUIP  yi/D  = 6,  r/R  = 1, 

0 = 135*;  OBSERVER  IS*’ia!  40*  TO  THE  JET  AXIS. 


FIG.  6.9  RELATIVE  COMIRIBUIION  TO  THE  SHEAR  NOISE  SPECTRUM  FROM  A 
UNIT  VOLUME  OF  JET.  SOURCE  POINTS  ARE  LOCATED  ALONG  THE 
LINE  r/R  =1,  0=0°. 
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FIG.  6.11  RELATIVE  CONTRIBUTION  TO  THE  SHEAR  NOISE  SPECTRUM  FROM  A 
UNIT  VOLUME  OF  JET.  SOURCE  POINTS  ARE  LOCATED  OVER  A 
CROSS-SECTION  AT  y^/D  = 5 • 
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FIG.  6.12  CONTOURS  OF  EQUAL  RELATIVE  CONTRIBUriON  FROM  UNIT  VOLUME 
TO  THE  OVERALL  SHEAR  NOISE  FOR  A SLICE  OF  JET  AT  yi/D  = 3 
NORMALIZED  BY  THE  CONTRIBUriON  FROM  y^/D  = 3,  r/R  « 1, 

0 = 0“j  OBSERVER  IS  AT  40*  TO  THE  JET  AXIS. 
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FIG.  6.13  CONTOURS  OF  EQUAL  RELATIVE  CONTRIBUTION  PER  UNIT  VOLUME  TO 
THE  OVERALL  SELF  NOISE  FOR  A SLICE  OF  JET  AT  yi/D  = 3 
NORMALIZED  BY  THE  CONTRIBUTION  FROM  yi/D  = 3j  r/R  = 1, 

0 = 0";  OBSERVER  IS  AT  40“  TO  THE  JET  AXIS. 
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FIG.  6.l4  NOMDIAL  SOURCE  PATTERNS  (ux^x^  ~ SHEAR;  (ux^)^  ~ SELF) 
AT  yi/D  “ 3.  THE  MAXIMA  ARE  FOUND  IN  THE  NEIGHBOURHOOD 
OP  r/R  = 1,  0 = 0*.  THE  VELOCZTY  COMPONENTS  IN  THE  x 
DIRBCTION  (40“  TO  THE  JET  AXIS)  HAVE  BEEN  MEASURED. 
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FK.  6.16  COlirOURS  OF  EQUAL  REUfflIVE  CONTRIBUTION  PER  UNIT  VOLUME  TO 
THE  OVIKALL  SHEAR  NOISE  FOR  A SLICE  OF  JEP  AT  yi/D  » 6 
NORMALIZED  BY  THE  CONTRIBUTION  FR(»1  y^/D  = 6,  r/R  - 1, 

<b  = 0*;  OBSERVER  IS  AT  40*  TO  THE  JET  AXIS. 
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FIG.  6.17  COOTOURS  OF  EQUAL  RELATIVE  COKERIBUT  ION  PER  UNCT  VOLUME  TO 
THE  OVERALL  SELF  NOISE  FOR  A SLICE  GF  JET  AT  yi/D  = 6, 

= 1,  0 = 0";  OBSERVER  IS  AT  40“  TO  THE  JET  AXIS. 
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FIG.  6.18  NCamiAL  SOURCE  PATTERNS  ~ SHEAR;  (ux^)^  ~ SELF)  AT 

jn/D  - 6.  THE  MAXIMA  ARE  FOUND  IN  THE  NEIGHBOURHOOD  OF 
r/E  » 1,  0 = 0*.  THE  VELOCin  COMPONENTS  IN  THE  x DIREC- 
TION (40“  TO  THE  JET  AXIS)  HAVE  BEEN  MEASURED. 
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FIG.  6.19  CONTRIBUTIONS  TO  THE  FAR-FIELD  SELF  AND  SHEAR  NOISE 

SPECTRA  FROM  A SLICE  OF  JET  AT  VARIOUS  AXIAL  POSITIONS 


FIG.  6.20  RELATIVE  OVERALL  SELF  AND  SHEAR  NOISE  SOURCE  STRENGTHS  PER 
UNIT  LENGTH  AS  A FUNCTION  OF  AXIAL  POSmON;  COMPARISON  OF 
EXPERIMENT  AND  A HYBRID  MODEL  PROPOSED  IN  APPENDIX  A.  THE 
DATA  HAS  BEEN  NORMALIZED  BY  THE  RESPECTIVE  VALUES  AT  yi/D  = 3. 
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FIG.  6.21  ACOUSTIC  SOURCE  STRENGTH  PER  UNIT  LENGTH  AS  A FUNCTION  OF 
FREQUENCY  AND  AXIAL  SOURCE  POSITION.  SOURCE  STRENGTH  DIS- 
TRIBUT  IONS  INFERRED  FROM  A POLAR  CORRELATION  SOURCE  LOCATION 
TECHNIQUE  (FISHER  ET  AL,  REF . 75)  ARE  SHOWN  FOR  QUALITATIVE 
COMPARISON. 
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FIG.  6.22  COMPARISON  OF  THE  FAR  FIELD  SHEAR  NOISE  SPECTRUM  COMPUTED 
FROM  CROSS_SPECTRAL  DENSITY  MEASUREMENTS  WITH  THE  ONE 
EXTRACTED  FROM  FAR  FIELD  JET  NOISE  DATA. 
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FIG.  6.24  TEST  OF  THE  SIMILARITY  OF  THE  PEAK  NORMALIZED  SELF  AND  SHEAR 
NOISE  SPECTRA  AFTER  A FREQUENCY  SHIFT:  BOTH  SHEAR  NOISE 
SPECTRA  HAVE  BEEN  SHIFTED  IN  FREQUENCY  BY  A FACTOR  OF  1.6. 


FIG.  6.25  OVERALL  FAR  FIELD  JET  NOISE  SPECTRUM  AT  40“  TO  THE  JET  AXIS: 

COMPARISON  OF  DIRECT  MEASUREMENT  WITH  PREDICTION  BASED  ON 
CROSS- SPECTRAL  DENSITY  MEASUREMENTS. 
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APPEMDIX  A 


A SIMPLE  JET  NOISE  MODEL 


Some  properties  of  jet  noise  are  cco^uted  from  a simple  turbulence- 
plus-mean  flow  model  which  is  largely  equivalent  to  that  of  Ribner  (Ref.  10) . 
The  analytical  approach  follows  along  the  lines  of  Pao  and  Lowson  (Ref.  78). 
By  transformation  to  the  wave  nxariber- frequency  space,  some  of  the  rather 
complex  coordinate  transformations  (see,  for  exaiiq>le,  Ribner  UTIA  Report  66, 
Ref.  6)  are  exchanged  for  more  straightforward  Fourier  transforms. 


The  turbiilence  model  that  will  be  used  to  describe  the  two  point 
velocity  correlations  in  the  jet  satisfies  the  criterion  of  homogeneous 
isotropic  turbulence  (Batchelor,  Ref.  79)  convected  by  the  local  jet  flow. 
The  cross-correlation  of  the  velocity  ui(x,t)  and  Uj(jj+^  t-r)  defined  in  a 
reference  frame  moving  at  the  convection  velocity  Up  is  modelled  as: 


Ui(^,t)Uj(y+^,  t-x)  = u^^(y)  |6^j(l-a^  ll^l^)  + 


-a2(  ||  |2  + a%^2) 
(Al) 


The  parameters  a,  a,  and  Ug  are  weak  functions  of  position  and  determine  the 
spatial  and  temporal  scales  of  the  correlation.  For  convenience  these  para- 
meters are  treated  as  being  constant  within  a correlation  volume.  The 
tvirbulence  model  is  identical  to  the  one  discussed  by  Pao  and  Lowson  (Ref. 

78) ; extensive  use  of  such  a model  was  also  made  in  earlier  work  by  Lilley 
(Ref.  5)  and  Ribner  (Ref.  10).  The  Gaussian  was  chosen  for  its  well  behaved 
Foxirier  transform;  even  after  repeated  differentiation  one  does  not  encounter 
convergence  problems.  Mathematical  expediency  has  dictated  the  choice  of  the 
model  correlation;  consequently,  the  tiirbulence,  which  is  neither  isotropic 
nor  homogeneous,  is  not  described  with  great  accuracy. 


The  model  cross- correlations  and  cross- spectral  densities  relevsuit 
to  Ribner 's  self  and  shear  noise  formalism  take  the  form: 


Ux\(y,^,T)  = u^‘^(^)  - a‘^(62Cose- lysine) 


■a^(  1 1 


(A2) 


U ^u'^(y>l>T)  = 2(u  U’)^  + (u  ^)^  (A3) 

X A XX  X 

Ux  Is  the  velocity  in  the  x direction  (the  observer  is  located  at  2^ , | the 
separation  between  field  points  in  the  turbulent  flow,  and  t is  the  tiS 
delay.  As  the  sovirce  term  of  the  jet  noise  is  proportional  to  32/^t2  (u*  + Ux)^> 
only  2(ux'Jx)2  in  Eq.  A3  will  contribute  to  the  jet  noise  correlation,  ( v^)2 
being  invariant  with  respect  to  t. 


A-1 


The  cross-spectra  in  wave-nuniber-frequency  space  are  given  by  the 
respective  four  dimensional  Foxirier  transforms  as; 

1 fu  |2  ^ 

(z)  ^ (k^cose  - iLsine)  4a  ^ a U 

„ (y,k,»)  = -V%  ^ 2 _S e 

l6a"au  r lla^ 


c 

2,  , ,2 


[\(y)  T 

<P  2 2^^’^’“) 1* 


u u 

X X 


32a  Q!U 


(a4) 

2 " “T  r 1^ 

r ^ ■ k^sine)  N,^2'j  8a  ^ 


u 
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k and  u are  the  wave  nurher  vector  and  the  radian  frequency  respectively. 
The  mean  velocity  profile  is  modelled  as: 


U 


o = 


2/  2 ^ 2a 

•c  (yg  + 73) 


(A6) 


with  c and  Uj_  slowly  varying  functions  of  axial  position  y]_.  The  mean  velocity 
ccn^jonent  in  the  direction  of  an  observer  at  x is: 


Ux  = Uo^ose 


(A7) 


The  Gaxissian  is  a reasonable  description  of  the  mean  velocity  profile  in  the 
transition  region  and  the  fully  developed  jet;  however,  it  fails  to  model  the 
flow  near  the  jet  nozzle  (mixing  region)  with  sufficient  acciiracy.  A better 
model  of  the  mean  velocity  profile  in  the  mixing  region  would  have  been  (Ref. 

80); 


00 
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(a8) 


Lilley  (Ref.  5)  suggests  A » 300  and  Jq  ^ *03 • Rihner  (Ref.  10)  models  the 
two-point  correlation  as 


U 


2/r2 

-CTlTTg  /L 
e 


(A9) 


rather  than  modelling  U itself. 

One  miist  keep  in  mind  that  the  correlation  fimctions  and  the  mean 
velocity  profile  are  but  crude  estimates  of  the  ones  that  would  be  measured 
in  a real  jet  flow.  The  spatial  and  ten5)oral  behaviour  of  the  model  fxmctions 
is  hopefully  illustrative  of  the  general  nature  of  the  'real'  functions. 


Jet  Noise  Estimate 


The  sheeir  noise  contribution  to  the  total  far  field  acoustic  pressure 
according  to  the  theory  of  Ref.  10  is: 


= "h,,  J ^ - V I*  - y i)<^  (aio) 


2vcflxl' 


The  shear  noise  auto- correlation  is  evaluated  from  a double  voluae  Integral: 


• ^ -Cl?-  ?i  l>  ^ "K(5e>‘  - - C li  - Sz 

(All) 

Ribner  (Ref.  10)  has  argued  that  cross- terms  such  as  [3^/St2(ibc)  ^/dt^(x3x^)  1 
which  are  found  in  a more  rigorous  description,  will  be  wman  or  identically 
zero.  After  transformation  of  the  variables  zi  = y - l/2  n,  Z2  = y + 1/2  ti, 
and  application  of  the  identity  ^ 


^ \ \ - TTT 


ht‘-  " 


hr 


u u' 

X X 


(A12) 


the  auto-correlation  takes  the  form: 

2 


o * 


2-,  y xv  2^  |x|c  / - - 


The  integration  over  evaluates  the  contrlbtitlon  to  Rsh(2^  t)  from  a correlation 
volume  at  It  is  convenient  to  replace  the  cross- correlation  by  the  equivalent 
Fourier  Integral  representation.  An  observer  moving  at  the  convection  velocity 
would  measvtre: 


^ ■/  V 
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-iUT 


dkdu) 


(A14) 


A-3 


Substitution  of  for  | changes  the  moving  frame  correlations  to  fixed 

frame  correlations; 


Vi  = j fu  u ’ 

'j  XX 


ik-(l-UT)  -i«T 


7.  rc 


dk  du 
-r 


(A15) 


It  follows  that 


C. 


o 'r' 


^-i(<o  - k-U^)  ^-i(k  + (u,  - k.U^)  f (1,32  + 


dy  dxi)  dk  dt) 

-T  ▼ T 


(A16) 


The  grouping  of  the  veiriables  in  the  exponentials  sv^ests  a change  of  variables: 


q = k + (u  - k'U  ) — 

^ 7 ^ TC^  C^  |X| 

(1 . (» - 


(A17) 


The  Jacobian  of  this  transformation  is  unity.  The  new  variables  will  sin?)lify 
the  integration  over  the  wave  number- frequency  space;  convection  effects  are 
also  accounted  for*  [viz  (u  - .l^'Uc)]*  Refraction  of  sound  has  been  siq)pressed 
by  replacing  the  pertiirbed  density  p within  the  jet  by  the  ambient  value  Pq 
(Ref.  24) . Integrating  over  jn  and  q-j^  results  in: 
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(AI8) 


where  C = [(1  - M^cose)  + ' is  the  convection  factor  first  found  by 

Ribner  (Ref.  5)  and  Ffowcs-Williams  (Ref.  11);  Mq  is  the  convection  Mach  nuiifcer, 
CqMc  = Uq  is  typically  half  the  jet  velocity  (Ref.  82) . After  integration 
over  q2  and  qo,  one  finds  that  the  shear  noise  auto- correlation  is  the  Fourier 
transform  of  the  spectral  contributions  from  eJ.1  correlation  volumes  in  the  jet 
flow; 
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(AI9) 


For  the  shear  flow  the  convection  factor  cones  out  to  he: 


Cgjj  = (1  - M^cose)^  + a\^  1 - ) 


(A20) 


There  is  virtually  no  difference  between  C and  Cgjj  at  low  convection  Mach 
numbers . 

The  volume  integral  of  ai i spectral  contributions  (deleting  the 
Fovtrier  transformation  in  Eq.  AI9)  defines  the  power  spectrxan  of  the  shear 
noise.  The  spectral  contribution  from  a unit  volxmae  is:  2 
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The  peak  frequency  of  dcpgjj  is 
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(A22) 


since  .5[a“^c^sin‘^e  • cos^e/(2a^  + c^)  ] is  much  smaller  than  (1  + cos"0) . 

The  mean  square  pressure  p „ = R (x,0)  can  be  written  in  the 


form: 


■SH  SH' 
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where  a = c^/a^;  a < 1,  as  the  correlation  length  (a  l/a)  cannot  be  greater 
than  the  Jet  diameter  (a  l/c) . QgH  interpreted  as  the  shear  noise 


contribution  per  unit  volvnne  to  the  overall  mean  square  pressvire.  The  inte- 
gral of  Osh  over  the  local  cross-section  of  the  jet  would  correspond  to  the 
shear  noise  contribution  to  the  source  strength  per  unit  length.  The  basic 
directivity  of  Qsh  is  proportional  to: 


2 ^ 
(l  + cos  e)  = cos* 


OV.- 


(A24) 


The  basic  directivity  is  modified  by  the  convection  factor  resulting  in 
a downstream  bias  and  refraction  which  is  not  included  herein.  Ribner  (Ref. 
10)  as  well  as  Pao  and  Lowson  (Ref.  78)  have  obtained  the  same  basic  direc- 
tivity for  the  shear  noise  term.  (Pao  and  Lowson  defined  the  shear  noise 
source  term  to  be 


5u 

SHEAR  NOISE  a sf:  '""5) 

which  appears  to  be  eqioivsilent  in  an  integral  sense  to  Ribner' s definition 
used  herein.) 

A similar  analysis  may  be  performed  for  the  self  noise  given  by: 

^O  '-7' 

The  mean  square  self  noise  pressure  is: 


(A27) 


with 
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The  spectral  distribution  is  proportional  to 
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jet  is; 


The  ratio  of  the  self  and  shear  noise  radiated  hy  a iinit  volume  of 


^ ^ (2  + ^ °SH  * 1 

U-^(y)  ?”  cos^0(l  + cos^e) 
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(A30) 


The  ratio  is  governed  by  the  spatisd.  scales  of  the  mean  flow  and  the  turbxilence 
(ff),  the  turbulence  intensity  distribution  and  the  basic  directivity  of  the 
shear  noise  (the  basic  self  noise  is  omnidirectional) . Convection  effects 
CsH^/c5  tend  to  cancel. 


The  self  and  shesur  noi,se  spectral  contributions  peah  at  different 
frequencies.  Nossier  and  Ribner(Ref.  37)  have  found  that  the  peak  self  and 
sheax  noise  frequencies  extracted  from  fax-field  data  differ  by  a factor  of 
about  2.  The  Gaussian  correlation  model  of  the  present  aneilysis  predicts  a 
peak  frequency  ratio  of  McCartney  (Ref.  77)  has  discussed  the  dependence 
of  the  frequency  shift  on  the  choice  of  the  tenporal  form  of  the  correlation 
and  cited  several  exanples.  The  model  self  noise  spectrum  can  be  modified  to 
exhibit  the  observed  frequency  shift:  The  argument  of  the  exponentieil  of  the 
self  noise  spectrum  c2/8a2a2Uc2  is  multiplied  by  a factor  of  l/2.  Hence  the 
peak  frequency,  defined  by  the  sensible  solution  of 
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becomes 
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Integration  over  fi  resiilts  in  a factor  of  25/2, 

thus 
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The  jet  half  width  is  given  by  th?  eii?)irical  relation  bW2  “ .0848y^  (Ref.  8l) ; 
thus  c = 9*62/yx*  The  turbulence  scale  is  estimated  to  ne  .L,=  .13yi  (Ref.  82), 
and  as  a = »/r/I<  the  square  of  the  ratio  of  the  length  scales  becomes  o-  » .52. 
cr  is  somewhat  sensitive  to  the  choice  of  X and  bi/2;  therefore,  the  veC-ue  of 
0-  = .45  used  by  Ribner  (Ref.  9)  hab  been  3iosen.  This  yields 
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a value  coo^atible  with  the  work  of  Nossier  and  Ribner  (Ref.  37)  • 


strength  Distribution  of  Noise  Sources  Along  a Jet 


The  functions  Qgp  and  QsH  describe  the  contribution  of  a unit 
volume  to  the  overall  self  and  shear  noise  respectively.  It  is  often 
assumed  that  the  tenporal  and  spatial  scales  of  the  turbulence  do  not 
change  over  a slice  of  Jet.  If  the  behaviour  of  ui^  Ul^  and  (ui2j2  gxe 
known,  one  may  estimate  the  contribution  of  a slic*''  of  Jet  to  the  overall 
Jet  noise. 

~~2 

The  indlvldvial  contribvttions  to  p from  \mit  volume  given  by  Eqs. 
an<i  A27  may  be  written  as; 
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Similarly, 
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Uq  is  a typical  frequency  that  is  assumed  to  scale  as  Ui(y]_)/y3_ 
(Ref.  1);  the  correlation  length  ff  varies  as  yq.  Thus  the  contributions  of  a 
slice  of  Jet  scale  as; 
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The  area  integrals  have  been  evaluated  from  measured  tiirbulence  intensities 
Vuj_^/Ui  and  mean  velocities  of  the  UTIAS  3/4"  model  Jet  (see  Fig.  4.6) . The 
functions  PsH(yi/D)/PsH(3)  and  PsF(yi/D)/PsF(3)  are  shown  in  Fig.  Al;  they 
are  proportion^  to  relative  self- shear  noise  contributions  to  the  far- field 
acoustic  power  from  a ;mit  length  of  Jet. 

Ribner  (Ref.  7)  and  others  (Refs.  5,  8)  used  similarity  and  scaling 
arguments  to  predict  that  the  source  strength  distribution  per  unit  length 
follows  a 1®’'^  (i.e.  constant  contribution  of  a slice  of  Jet  to  the  far 
field  acoustic  power)  in  the  mixing  region,  and  a y“^  law  in  the  fully  developed 
Jet.  The  present  model  is  a hybrid  combining  scaling  considerations  (uo  ot  U^/yij 
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L O yi}*and  certain  products  of  measvired  jet  flow  velocities,  as  described 
by  Eqs.  A36  and  A37>  The  peak  turbulence  intensity  is  not  constant  in  the 
mixing  region  reaching  a maximvmi  vedue  near  y^/D  = 5.  Hence  Ribner's  yi° 
law,  predicated  on  the  assuiqptlon  of  constant  tvirbulence  intensity,  is  not 
satisfied.  Downstream  of  the  potential  core  (yi/D  >4)  Psp  and  Psh  decrease 
rapidly,  and  do  approach  the  yi"^  law  predicted  for  the  fully  developed  jet. 

Hypothetical  Cross-Correlations 

The  procedvire  of  estimating  the  cross-correlations  Uif(>,t)p(x.t-T^ 
and  U3f^(y.t^p(x.t-TV  is  similar  to  the  one  employed  in  the  cooptation  of  the 
far  field  pressure.  The  analysis  evaluates  the  cross- correlation  functions 
that  would  be  measured,  if  the  flow  were  accurately  described  by  the  assumed 
correlations  of  the  model. 

The  shear  noise  correlation  is: 


= u (y,t)p(x,t  + t) 
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where  ^ is  the  separation  of  an  arbitrary  point  and  the  source  point  y. 

The  cross- correlation  Uxui  is  pon-zero  only  for  small  values  of  |||  (corre- 
lation volume  of  Ux) . The  variation  of  the  mean  velocity  Ux  over  the 
correlation  volume  will  distort  the  iptegrand.  For  the  case  of  homogeneous 
isotropic  turbulence  and  flow  parallel  to  the  jet  axis,  the  Integrand  will  be 

maximum  at  jr  + e,  where  e is  a slight  offset  toweirds  the  jet  centre  line. 

Upon  substitution  for.Ux~and  replacement  of  the  cross  Correlation  ^/br^( uxUx) 
with  the  equivalent  Fourier  integral,  R\j^(x,y-,T)  becomes: 
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After  the  change  of  variables: 
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one  observes  that 
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is  the  three  dimensional  Fourier  transform  of  e 


and 
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A lengthy  but  straightforward  integration  results  in: 
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where 
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The  'interaction'  of  the  mean  shear  and  the  two  point  turbulent 
velocity  correlation  has  generated  two  functions,  one  of  which  is  even  about 
T = 0,  the  other  odd.  In  effect  the  mean  shear  has  destroyed  the  homogeneity 
and  the  isotropy  of  the  correlation.  The  turbulence  structure  of  the  jet  flow 
is  not  homogeneous  and  isotropic,  and  one  would  ejqject  the  measured  cross- 
correlations of  the  shear  noise  source  term  with  the  far  field  sound  to  be 
somewhat  sensitive  to  the  position  of  the  source  point  in  the  jet. 

The  relative  asplitudes  of  the  even  and  odd  functions  predicted 
above  are  plotted  in  Fig.  A2.  The  apparent  position  of  the  source  has  shifted 
from  y to  (y^,  (l+a)"^  y2,  y3) . The  shift  is  not  very  significant  from  the 
point  of  view  of  the  percentage  change  in  arrival  time.  The  shift  may  became 
significant  if  axial  variations  in  the  mean  flow  (neglected  herein)  are 
considered. 

In  order  to  conpute  the  contribution  to  the  far  field  shear  noise 
frcm  a unit  volume  at  y,  one  must  double  differentiate  the  cross-correlation 
with  respect  to  t andwaduate  the  result  at  t = Cq"^  (in  the  frequency 

domain  the  cross- spectrum  must  be  multiplied  by  -(w-h-ijc)^) . The  apparent 
source  shift  isplies  that  the  resulting  function  is  evaluated  near,  but  not 
at,  its  maximvon.  Should  one  neglect  the  small  time  shift  one  would  overestimate 
the  shear  noise  contribution  (the  same  applies  to  the  estimate  of  the  far  field 
shear  noise  spectrum;  here  the  time  shift  takes  the  form  of  an  additional  phase 
angle) . The  error  can  be  considerable,  as  the  double  time  derivative  brings 
about  a narrower  peak. 


The  cross-correlation  Rujjp(x,y,T)  may  be  normalized  by  the  RtC  values 
of  Ux  and  p which  are  given  by: 


A-11 


T 


for  a = .45,  L = <Jv  a"^,  0 = 40“. 

yjEi;/L3  has  been  interpreted  (Ref.  39)  as  the  n;in4)er  of  uncorre- 
lated volume  soitrces  assianed  to  be  of  equal  strength;  the  correlation 
coefficient,  therefore,  is  inversely p roportional  to  the  square  root  of  the 
number  of  uncorrelated  sources  in  the  jet. 

The  self  noise  correlation  is; 
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does  not  suffer  from  the  shear  flow  interaction.  After  some  algebra; 
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The  normalized  self  noise  correlation  is  given  by; 
1 1 

when  cr  = .45,  L = ^y^■  a”^,  0 = 40®. 


Summary 


A sin5)le  turbulence-plus-mean  flow  model  has  been  used  to  illustrate 
what  is  hoped  to  be  the  general  nature  of  the  measured  functions  Ru^^p  and  Ru^^p* 

The  model  predicts  that  the  peak  of  the  shear  noise  correlation  Ru^^p  will  occur 
at  a time  delay  different  from  the  travel  time  c^^  1*"^  !•  Measvired  Ruj^p  cross- 
correlations (see  Section  6.1)  indicate  that  the  predictions  considerably 
underestimate  the  shift  in  the  time  delay. 


A- 12 


RELATIVE  RELATIVE  SOURCE  STRENGTH 


-.1.0 


2 3 4 5 6 7 y,/D 

FIG.  A.l  RESiATIVE  SOURCE  STRENGTH  PER  UNIT  LENGTH  AS  A FUNCTION  OF 
AXIAL  POSITION.  THE  CURVES  HAVE  BEEN  NORMALIZED  BY  THEIR 
RESPECTIVE  VALUES  AT  yi/D  = 3- 
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APPENDIX  B 


NOISE  GENERATED  BY  TRANSDUCER-FLOW  INTERACTION 


Even  though  the  principal  velocity  measuring  system  used  dxiring 
the  present  investigation  was  the  Laser  Doppler  Velocimeter,  sene  thought 
Was  given  to  the  problem  of  transducer- flow  interaction,  as  the  hot-wire 
anemcneter  will  continue  to  be  the  most  common  velocity  sensor  used  in  fluid 
flow  research.  The  instrument  is  aarailable  ’off  the  shelf  and  does  not 
entail  the  extensive  ancillary  system  development  reqvilred  for  the  Laser 
Doppler  Velocimeter.  Fvirthermore , the  volume  of  the  active  sensor  is  very 
smal .1  conpared  to  the  measuring  volume  of  a Laser  Doppler  Velocimeter, 
making  the  hot-wire  anemometer  ideally  suited  to  survey  scale  model  flows. 

A solid  object,  such  as  a hot  wire  probe  and  its  sx^jports,  radiates 
sound  when  placed  in  a txirbulent  flow.  Curie  has  provided  an  analytical 
framework  that  allows  one  to  predict  the  soxind  generated  by  such  flow-surface 
interactions  (Ref.  k2) . The  theory  identifies  fluctuating  forces  (lift  for 
exanple)  as  sound  generators.  This  is  confirmed  by  direct  correlation  of  the 
force  with  the  radiated  sound  as  was  done  by  Clark  and  Ribner  (Ref.  43)  • 

Cxirle  shows  that  dipole  radiation  due  to  flow-surface  interaction  can  be  of 
ccoparable  intensity  with  the  quadripole  radiation  of  a txirbulent  flow,  if  the 
Mach  number  is  small  enough. 

A hot  wire  probe  when  properly  placed  in  a Jet  flow  does  not  produce 
a measurable  increase  in  the  radiated  noise  (it  is  unlikely  that  a 1^  increase 
in  sound  pressure  l.ivel  can  be  detected  with  any  confidence) . In  Appendix  A 
it  is  shown  that  the  nomalized  crosi^ correlation  coefficients  of  the  velocity- 
pressure  correlations  are  of  order  if  there  are  N uncorrelated  radiators  of 
equal  strength  Each  soiuce  contributes  an  {onomit  to  the  mean  square 
pressure.  The  mg^  square  pressure  generated  by  the  probe- flow  interytion 
assumed  to  be  e = 0(l).  The  new  overall  RJC  sound  pressure  isV^N+^^S^ 

and  as  N » € there  is  little  change  in  the  far  field  pressure.  The  correlation 
function  of  the  measured  velocity  amd  the  Jet  noise  is  now  of  order 
as  the  sound  generated  by  the  probe  is  well  correlated  with  the  measured  velocity, 
the  turbulent  flow  being  responsible  for  the  fluctuating  forces  on  the  probe. 

Since  6 = 0(1)  the  error  due  to  probe  noise,  which  is  an  artifact  of  the  measure- 
ment, can  seriously  'contaminate'  or  obscure  the  Jet  noise-jet  flow  cross- 
correlation . 

The  interaction  of  a sound  wave  with  a solid  object  generates  a 
diffracted  wave  by  virtue  of  the  reaction  of  the  surface  to  the  forces  exerted 
on  it  by  the  incident  wave.  If  the  body  is  at  rest,  the  reaction  is  such  as 
to  maintain  zero  velocity  normal  to  the  surface.  The  velocity  distribution 
necessary  to  cancel  the  incident  normal  velocity  can  be  generated  by  a collection 
of  acoustic  sources  located  inside  or  on  the  sxirface  of  the  body  (Morse,  Ref.  83). 

Similar  boundary  conditions  are  encountered  when  a body  is  immersed 
in  a turbulent  flow.  The  forces  exerted  on  the  body.  Just  as  the  stresses  in 
the  turbulent  flow,  act  as  sound  sources.  In  the  presence  of  solid  boundaries 
the  solution  of  the  wave  equation. 
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is  given  by  (Ref.  4?): 
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where  are  the  direction  cosines  of  the  outwaurd  drawn  normal  from  the  fluid. 
Equation  B1  is  an  approximate  description  of  sound  generation  by  a turbiolent 
flow.  Curie  (Ref.  42)  has  applied  the  divergence  theorem  two  times  to  the 
volume  integral  and  found  that  for  the  special  case  of  a stationary  surface 
the  radiated  sound  pressure  is; 

1 r ^2  6(t'  - t + c'^lx  - ;(|) 

1 r 5k  6(t'  - t + c"^lx  - y|) 

■ ^ 
s 

where  Fj[^  = /ibijP  is  the  force  per  unit  area  exerted  on  the  fliaid  by  the  solid 
boundaries.  In  the  far-field:  ' 

P(x,t)  ^ r£-v/(y,t')  6(t'  .t  + c;^Ix-yl)  dy  dt' 

+ F (y,t’)  6(t'  - t + c‘^|x  - j^l)  dy  dt'  (b4) 

4ir:  IxlJ  o -r  ^ ^ 

o •*»  * 

The  Reynolds  stresses  pv^vj  and  the  forces  have  been  replaced 
by  the  momentum  flux  Po^x^  a^id  a force  Fx  the  direction  of  the  observer 
at  If  there  is  no  fluctuating  force  in  the  direction  of  the  observer,  then 
the  observer  will  hear  no  sound  due  to  flow-surface  interaction.  This  can  be 
illustrated  by  the  sinqple  exan^jle  of  an  oscillating  sphere  (a  dipole)  set  in 
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motion  by  a force  in  the  horizontal  plane.  No  soxmd  is  radiated  in  the 
direction  normal  to  the  dipole  axis  which  coincides  with  the  force  vector. 

The  forces  exerted  on  the  body  in  a turbulent  flow  are  of  aero- 
dynamic origin.  Rackl  (Ref.  84)  reasons  that  the  force  can  be  modelled  as: 
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where  the  pure  number  K is  a function  of  the  body  geometry  and  S is  the 
surface  area  exposed  to  the  flow.  The  far  field  acoustic  pressure  beccmes: 


P(x>t) 
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provided  the  scale  of  the  body  is  small  compared  with  the  scale  of  variation 
of  Vx. 

The  cross- correlation  of  Jet  flow  with  the  radiated  Jet  noise  gives 
rise  to  two  types  of  correlations:  one  associated  with  the  sound  generated 
by  the  turbulent  flow  alone,  the  other  generated  by  flow-probe  interaction. 

The  relative  magnitudes  of  the  two  correlations  can  be  estimated  with  the  aid 
of  a simple  model.  As  was  the  case  with  the  self  and  shear  noise  model,  v 
separates  into  a mean  velocity  IL.  and  a turbulent  velocity  Uvfv.t) . The 
cross-correlations  Ru^p  = Ux(y,tTp(x,t+T)  and  Ru^^p  = Ux^(y,t)p(x,t+T)  are 

to  be  estimated.  The  measured  B\:ijj2p  correlations  are  too  poorly  defined  to 

be  analyzed  with  any  confidence;  therefore,  only  the  R„  p correlation  will  be 
examined . 

• 

The  velocity  is  well  correlated  over  the  prove  sensor  and  the  portion 
of  the  probe  supports  that  are  immersed  in  the  correlation  volume;  hence,  the 
cross-correlation  is  given  by: 
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^eff  ® measure  of  the  surface  area  of  the  probe  in  the  correlation  volume. 

The  analysis  in  Appendix  A shows  that  the  convection  effects  are 
at  low  Mach  numbers;  we  ^all  neglect  it  here.  Furthermore,  the  major 
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effect  of  the  mean  shear  is  to  alter  the  directivity  factor;  hence  it  is 
assvimed  to  be  constant  over  a correlation  volxme.  Th\is  the  rather  siiqple 
model  correlation: 

^ = \ ^ 3 (Be 


is  chosen  for  mathematical  convenience.  Upon  substitution  in  Eq.  B7,  the 
smal 1 Mach  nunfcer  approximation  of  the  cross- correlation  takes  the  form:. 
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RujjP  is  the  sum  of  an  even  fxmction  and  an  odd  fvinction.  The  latter 
is  the  'contamination'  caused  by  the  probe-flow  interference  auid  must  be  small 
if  the  jet  flow- jet  noise  correlation  is  to  be  estimated  with  reasonable 
accuracy.  The  ratio  of  the  maximum  valxoes  of  the  two  correlations  will  g;ive 
an  indication  of  the  iu^portance  of  the  probe  noise.  The  even  function: 
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has  a maximum  at  t = + •^(j2a(l-*a)V2u^.  The  ratio  of  the  maximum 


a^Sg^VlT? 
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should  be  large,  if  the  probe  noise  is  negligible.  This  criterion  may  be  met 
if  the  effective  surface  area  is  amaii  and  the  Mach  nunher  relatively  high. 

Let  L be  a typical  correlation  length;  it  follows  that  a^  = tt/l^ 
and  Seff  = ^tflpL,  where  Rp  is  the  radius  of  the  probe  support  needle  (the  hot 
wire  sensor  is  much  smaller  than  the  probe  support) . The  correlation  length 
scales  with  the  axial  position  in  the  Jet:  L = .13  B (yi/B),  D is  the  diameter 
of  the  Jet  nozzle.  The  Reynolds  nusber  of  a cylindrical  probe  support  of  radius 


.01"  in  a Mach  .3  air  Jet  is  appraxinately  eqiial  to  3^00,  and  the  drag 
coefficient  is  near  unity  (Schlichting,  Ref.  81) . The  parameter  K is  set 
to  be  equal  to  the  drag  coefficient,  and  a is  chosen  to  be  .55  (Ref.  10) . 

For  a Jet  diameter  of  3/4”  and  Uq  = l/2  U , the  ratio  of  the  mAarlnm.  is; 

.16  ~ (BI3) 

The  model  predicts  that  the  probe  noise  makes  a significemt 
contribution  to  the  correlation  function,  when  the  probe  is  placed  in  the 
principal  sound  producing  region  of  a Jet  flow,  nanely  the  first  15  Jet 
diameters  downstream  of  the  nozzle.  Such  a prediction  is  conqpatible  with  the 
marked  differences  of  the  shear  noise  correlations  measured  with  the  Laser 
Doppler  Velocimeter  and  the  ones  measured  with  the  hot  wire  anemone  ter. 


APPENDIX  C 


CROSS-SPECTRAL  DENSITY  ESIIMAIES 
COMPUTED  BY  A FAST  FOURIER  TRANSFORM  ANALYZER 


The  fast  Fovirier  transform  (FFT)  analyzer  exploits  the  close 
relationship  between  the  Fovirier  series  coefficients  of  a periodic  wave 
"S(t)  and  the  Fourier  transform  of  the  non-periodic  version  of  the  same  wave 
shape,  namely  that  the  Fourier  series  coefficients  are  unifonnly  spaced 
sanples  of  the  Fovirier  transform  of  the  non-period  wave  (Harils,  Ref.  85). 
For  the  pvirpose  of  machine  confutation  only  a finite  number  of  points  at 
which  s(t)  is  known  can  be  used,  and  similarly  there  will  be  a finite  nianber 
of  frequencies  at  which  the  Fovirier  transform  of  the  truncated  sequence  of 
S(w)  can  be  conputed.  The  fast  Fovirier  analyzer  estimates  the  Fourier 
transform  of  an  arbitrary  fvinctlon  f(t) , subject  to  the  constraints,  imposed 
above. 

The  effect  of  the  finite  duration,  dviring  which  f(t)  is  being 
analyzed,  is  equivalent  to  passing  the  signal  f(t)  through  a time  window 
D(t): 


D(t) 


r T 

0,  t < - •^ 

■ Seme  suitable  function 
T 

0,  t >1 


(Cl) 


D(t)  is  chosen  to  minimize  the  errors  inherent  in  the  process  of  power  spectreLL 
density  estimation.  The  most  common  errors  eure  spectral  broadening  (finite 
analysis  time)  and  the  generation  of  side  lobes  (dlscontlnviitles  at  the 
beginning  and  end  of  the  window) . The  behaviour  of  time  window  has  been 
discussed  in  some  detail  by  Harris  (Ref.  86) . 

The  Fovirier  transform  of  the  function  D(t)f(t)  is 


^(u)  D(t)f(t)e^"*^  dt  (C2) 

-00 

Consider  the  product  of  the  Fourier  transform  estimates  of  f(t)  and  g(t): 


=-^5  rD(tjf(t)e^"V(f)g»(t’)e'^"*''  dt  df  (C3) 

Utt  J 

The  window  fvinctlon  has  a finite  duration  and  possesses  a Fovirier  transform 

JCr(«); 


c-1 


and 


D(t)  = J'  «5l(«)e  du 


4 /'  , r -i(«  -fcL)t  i(«L-«)T 

^ (u)^(u)  = ^(u^)^(u^)f(t)gx-(t  + T)e  ° ^ e ^ dt  dw,  d«  dt' 

4Tr 

(C5) 

A A 

Each  realizat^n  of  3^  * can  be  treated  as  a mentoer  of  an  infinite  set  of 

realizations  “oi  S*)  if  the  process  is  ergodic.  The  ensemble  average  is 
defined  to  be; 


<'i?i'«J?'i>  = qJffCw)  = <fi(t)g?(t  + t)> 

4ir  o 


-i(w^-u^)t  i(«^-tt)T 


dt  diij^  du^  dt' 


X 

<fi(t)gj(t  + t)>  = R^g(T)  = ^ J f(t)e*(t  + T)dt  (C7) 

-T 

(Blackman,  Tukey,  Ref.  87),  and  the  correlation  fvinction  is  independent  of  t. 
The  integration  with  respect  to  time  t will  result  in  the  delta  function 
2v  6(«q  - «^) ; thus 

^fg^“^  "f  l^“i)l^  '’’fgK  ■ “)  <401  (C8) 

where  <?^g(u)  is  the  cross- spectrsJ.  density  of  f(t)  and  g(t) . The  convolution 
integrsd  accounts  for  the  absence  of  all  |t|  > T,  brought  about  by  the  finite 
window  duration. 

5"  ^ analyzer  computes  the  average  value  of  N realizations  of 

w i^(«) . Each  average  may  then  be  treated  as  an  estimate  of  ^ (w). 
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APPENDIX  D 


SOME  PROPERTIES  OF  THE  LASER  DOPPLER  SIGNALS 


wave 


A laser  beam  travelling  in  the  direction  can  be  modelled  by  the 


= E^(xj^)e 


2^  2. 
x„  4x  ) 


i(2w^t. 


kxi) 


(Dl) 


Vq  is  the  frequency  of  the  laser  light  with  a wave  nimfeer  h.  Eq  and  a are 
slowly  vaxylng  functions  of  xi  and  account  for  the  slow  beam  divergence. 
When  two  beams  derived  frcm  a conmaon  soiirce  intersect  at  x = 0 with  the 
angle  20,  the  total  field  near  x = 0 is 

- a^  ( X-  ^ sin^04Xo^  cos^0+X3^ ) i ( 2 ttw  t-loc,  CO80) 
i^(x,t)  = 2E^  e ^ 2 3 , 


2 

•cosh(2a  x^X2cos0sin0+ikX28in0)  (D2) 

2 

The  intensity  of  the  light  is  proportionsQ.  to  IBj  | . The  intensity  distribution 

2 2 2 2 2 2 
o o ■2a  (x,  sin  0^x_  cos  0+x_  ) 

|E^r  = 2E/e  12  3- 

O 

^ (cosh(4a  x^X2Cos0sin0)  + cos(2kX2sin0) } (D3) 


is  the  sum  of  a smooth  function  known  as  the  Doppler  pedestal  and  an  inter- 
ference pattern  cos(2kx2sin0) . The  distance  between  two  maxima  in  the 
Interference  pattern  is  the  fringe  spacing 


r 2sin0 


(D**) 


The  Gaussian  intensity  profiles  of  the  laser  beams  confine  the  fringe  pattern 
near  the  point  of  intersection  x = 0 ususdly  referred  to  as  the  probe  volume 
(Fig.  Dl) . Surfaces  of  constant  fringe  intensity  are  defined  by 

2 2 2 2 2 

Xj^  sin  9 + Xg  cos  0 + x^  = constant  (D5) 

The  probe  volume  is  in  the  shape  of  an  ellipsoid.  The  interference  fringes 
are  paredlel  to  the  X1X3  plane  and  are  separated  by  a distance  df.  The  angle 
of  intersection  20  is  usually  in  the  range  of  3*  to  40**,  and  the  ellipse  in 
the  X1X3  plane  has  a large  eccentricity,  whereas  the  ellipse  in  the 


D-1 


d 


plane  is  nearly  circular.  For  sufficiently  large,  the  Doppler  pedestal 
meisks  the  fringe  pattern,  and  the  Doppler  signal  that  is  scattered  hy  particles 
can  no  longer  be  detected.  Signal  to  noise  considerations  limit  the  LDV 
processor  to  accept  only  Doppler  signals  that  exceed  a certain  threshold  level 
and  the  probe  volume  is  bounded  by  the  svirface  (Fig.  D2) 

x^^sin^0  + Xg^cos^^  + (d6) 

The  processor  further  requires  that  the  scattering  particle  cross 
at  least  No  fringes  (Nq  is  usually  8) . A particle  that  traverses  the  probe 
volume  neax  the  outer  boundary  may  not  intercept  a sxjfficient  number  of  fringes. 
The  confutation  procedure,  already  initiated,  cannot  be  ccmpleted.  The  per- 
centage of  valid  Doppler  bursts  in  the  total  nuntoer  of  signals  that  initiate 
confutation  is  defined  to  be  the  ^ validation.  The  largest  possible  value 
of  the  % validation  is  the  ratio  of  the  projected  areas  of  the  effective  probe 
volume  containing  at  least  8 fringes,  and  the  actual  probe  volume.  The  analysis 
applies  only  for  laminar  flow.  The  presence  of  turbxxlence  will  degrade  the  % 
validation. 

There  is  a total  of  N = 2R/df  fringes  in  the  probe  volume,  and 
therefore  the  effective  probe  volume  has  semi -major  axis 

a,  . (/  - N (D7) 

^ 2sin0  ' ' 


and  semi -minor  axis 

b^  = (1^  - (D8) 


The  area  ratio  is 


(D9) 


100  (1  - No^/N^)  is  the  largest  possible  % validation.  To  make  the 
best  use  of  the  scattering  particles,  the  % validation  should  be  near  the 
theoretical  maximum.  The  % validation  is  an  indirect  measure  of  the  data  rate. 
The  data  rate  fc  depends  on  the  number  of  particles  that  cross  the  prctoe 
volume  per  unit  time: 

n/  \ 

X-f)  (MO) 


If  a velocity  conponent  other  than  the  X2  direction  is  to  be 
measured,  the  probe  volvuae  must  be  rotated  about  the  x^  or  X3  axis.  As  the 
probe  volume  is  rotated,  the  scattering  cross-section  of  the  prcbe  volume 
changes.  If  the  scattering  cross-sections  of  the  probe  volume  and  the 


D-2 


1 


r 


effective  probe  volume  change  at  the  same  rate,  then  the  ^ validation  will 
be  unaffected.  The  data  rate  changes,  as  the  number  of  particles  that  cross 
the  effective  probe  volume  Is  a function  of  the  scattering  cross-section. 

Rotation  about  the  or  X2  axis  is  analogous  to  the  case  of  laminar 
flow  malting  an  angle  9 with  respect  to  the  normal  of  the  interference  fringes. 
Rotation  about  the  xi  axis  is  considered  first.  The  particles  that  intercept 
the  probe  volume  (Fig.  D3)  move  along  lines  defined  by: 

y = - (z  - z)tan0  (Dll) 

2 2 2 2 2 2 2 

and  intersect  the  ellipse  y cos  0 + z = R - x sin  0 = R (x)  at: 


A 2 2 2 2 2 

z tan‘^9  cos‘^0  ± TR  (x)(1  + tan  9 cos  0) 


'^2^  2. .1/2 
z tan  9 cos  01 


1 + tan  9 cos  0 


R(x)(1  + tan^9  cos^0)^'^^ 
tan9  COS0 


(D12) 


*2  - ^1  = 


p o 2 2 '‘2  2 2 l/2 

^2 ^ [R  (x)  1 + tan  9 cos  0 - z tan  P cos  0]  ' 


1 + tan  0 cos  0 


By  inspection: 


I K = cos20)^/2 


(D13) 


The  scattering  cross-section  of  the  probe  volume  is  an  ellipse  with  semi-minor 
axis: 


>>0 ' 5 


= R(1  + tan^e  cos^0)^'^^ 

COS0  ' 


(D14) 


Since  the  probe  volume  is  rotated  about  the  x^  axis,  the  semi-major  axis  is 
Invariant: 


^■0  ~ sin0 


(D15) 


The  scattering  cross-section  of  the  effective  probe  volume  is  bounded  by  the 
locus  of  all  points  such  that: 


*2  “ *1  “ ”o*^f 
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(D16) 


L. 


and 


i (1  + tan^e  coB%)^^  (1  + tan^e  cos^0)^/^ 

The  semi -minor  axis  is  defined  by  l/2  i(0)  and  the  semi -major  axis  has  length 
X,  \diere  jl{k)  = 0: 

1 / 2 2 2 

- -nr^  (1  + tan'^e  cos‘'0)J  (D17) 

2 2 

" coif  TT^  cos^0)^  (l  + tan^e  cos^0)^/^  (Dl8) 

The  area  ratio 


N 


= 1 - 


(1 


. 2 
tan  e 


COS^0) 


As  0 is  increased,  a critical  angle  0crit  reached  at  which  the 
area  ratio  is  zero.  At  this  particular  angle  (dead  angle)  and  beyond,  no 
valid  velocity  measurements  can  be  made  with  the  LDV  system; 


0 ..  = arctan 

crit 


(D19) 


A similar  auialysis  can  be  performed  for  the  case  of  rotation  about 
the  X3  axis  (Fig.  d4).  The  resxilts  are: 

No^  2 2 

area  ratio  = 1 - (l  + tan  0 sin  0) 

(D20) 

0crit  = arctan 

o 

In  many  systems  the  angle  0 is  of  the  order  of  a few  degrees,  and 
the  latter  configxiration  offers  a better  % validation,  whereas  the  former 
possesses  a sv^erior  data  rate  (Figs . , D6) . 

In  the  present  investigation  rotaticxi  about  the  X3  axis  was  chosen, 
even  though  this  particular  configuration  is  not  the  best  from  the  point  of 
view  of  maximum  achievable  data  rate.  The  system  does,  however,  make  use  of 
nearly  all  particles  that  cross  the  probe  volume  and  it  is  less  susceptible 
to  mean  velocity  changes  across  the  probe  volume. 
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NORMALIZED  DATA  RATE  ►.d  % VALIDATION 


. D.5  1o  VALIDAIION  AS  A FUNCTION  OF  ERQBE  VOLUME  ORIENTATION  WITH 
RESPECT  TO  THE  FLOW  DIRECTION. 


FIG.  D.6  data  rate  as  A FUNCTION  OF  PROBE  ORIENTATION  WITH  RESPECT 
TO  THE  FLOW  DIRECTION. 


APPENDIX  E 


ESTIMATE  OF  THE  POWER  SPECTRAL  DETSin  OF  THE  LDV  OUTPOT 


The  laser  Doppler  veloclmeter  generates  a sequence  of  ou^ut  voltages 
v(Tn)  proportional  to  the  speed  of  a particle  in  the  probe  volume  at  time 
Tn  - 5>  6 is  the  conputatlon  time  required  by  the  processor  (Fig.  El).  The 
last  realization  of  v(Tn)  is  stored  in  the  output  memory  of  the  processor  until 
a new  measurement  has  been  performed.  The  time  interval  between  successive 
measurements  v(Tn+i)  and  v(Tn)  is  AEn>  and  the  output  waveform  of  the  LDV  has 
the  form 

N 

^(t)  v(T^)[H(t  - T^)  - H(t  - (El) 

n=-N 

The  Fourier  transform  of  v(t)  is 
^ N 

^(«)  = ^ v(T^)[H(t  - T^)  - H(t  - ^1)16^*^  dfc  (E2) 

n*-N 


N 


V- 


itjr, 


iw 


h 


n e 
n 


n=-N 


sin 


(E3) 


u 


For  the  special  case  of  /!Tn  = define  a sazopllng  frequency  or  data 

rate  fc  = 


(sin  T?f/fc  ® ^ ®®n  be  interpreted  as  the  transfer  function  of  the 

particular  sampling  process.  The  asplltude  transfer  function  (sin  77f/fc.)/^/^c 
has  zeroes  at  f = nfc>  n = 1,  2,  ...  . The  frequency  window  is  a good  low  pass 
filter  for  frequencies  less  than  f fc/^  (FiS*  SS) . The  summation  is  an 
estimate  of  the  Fourier  transform  of  v(t) ; 


sin  77f/fj,  iTif/f^ 


N 


lun^  _ 
e ^ 


(e4) 


n*-N 


^^(«)  Jv(t)  e^*^  dt  (E5) 

D" (f)  is  a good  estimate  ofZf^f)  if  f < fc/2  and  if ^(f)  has  little 
spectral  content  for  frequencies  greater  than  fc/2.  If  this  is  not  the  cause, 
then  the  estimated  spectrum  will  contain  spurious  frequency  conponents,  called 
allauses. 
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As  a sixqple  excaqple  one  nay  consider  the  case  of  oscillatory  flow 
with  frequency  fQ.  The  output  of  the  LDV  processor  will  he; 

N 

W(t)  = cosWoTnt”^*  - ^n^  - H(t  - ] (e6) 

n=-N 

^ sin  77f/f  i7rf/fc  f 1 n ^ 

= ^o  “Wf—  ® J-  e (E7) 

^ y^=_YI  ^ 


sinirf/f  iT7f/fc  f i ^ n ^c  ^c  M 

n=0 


The  series  can  be  written  in  the  form; 


N f. 


^I( 


C nU 


+ e 


f-f 

2T7i 

c \n 


+ e 


T^( 


c \n 


For  N ->oo  the  series  simplifies  to  (Gradshteyn  and  Ryzhik,  Ref. 

. « 4 ^ 


(ElO) 


1 - e 


1 - e 


1 - e 


1 - e 


= l/27r  unless  f = ± f^  or  f = mf^  - ^o’  “ integer. 


At  these  special  frequencies  the  series  has  an  infinite  value,  and  the  usual 
interpretation  of  a delta  function  applies.  The  function  W(f)  has  maxima  at 
f = ± fjj,  the  frequency  of  oscillation,  and  at  f = m fg  ± fg.  The  latter  are 
the  aliases.  There  are  no  aliases  in  the  range  0 < f < fc/2  for  the  case  of 
fo<V2- 

In  order  to  avoid  aliasing  errors,  the  data  rate  of  the  LDV  must  be 
sufficiently  high  and  satisfy  the  criteria  outlined  above.  Turbulence  spectra 
of  Hot  Wire  Anemometer  signeds  indicate  that  fg  should  be  20  KHz  or  greater  if 
aliasing  errors  are  too  negligible.  Additional  filtering  of  the  signal  v(t) 
generally  does  not  reduce  errors,  as  the  aliases  may  be  distributed  over  the 
entire  span  0 < f < f^. 
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FIE.  E.2  SPECTRAL  WINDOW  IMPOSED  ON  THE  SPECTRUM  ESTIMATE  OF  THE 

OUTPUT  OF  THE  ILV  PROCESSOR,  AS  A CONSEQUENCE  OF  THE  DATA 
RATE  fg. 
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partially  successful.  The  major  difficulty  has  been  suspected  as  spurious 
•probe  noise'  generated  by  turbulence  - hot  wire  interaction.  Thus,  to  avoid 
this  problem,  the  traditional  hot  wire  anemometer  has  been  replaced  in  the  presen 
investigation  by  a non-intrusive  device:  a Laser  Doppler  Velocimeter.  Substanti 


UD  I JAN  73  1473  IIM('I  ACCTITTL'n 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fH7>*n  Dtim  Bnitttd) 


5ECUHITV  classification  OF  THIS  HAGEfHTiMi  /)•!•  KnIfrrJ) 


modifications  were  made  to  meet  tlic  constraints  i+nposcd  by  the  correlation 
experiment;  a major  feature  was  provision  to  measure  u^,  the  component  of 
turbulent  velocity  in  the  observer  direction  x.  Cross-correlations  and  cross- 
spectral  densities  of  the  jet  noise  at  40°  to  the  jet  axis  and  the  instantaneous 
turbulent  jet  flow  Uj^  Ux/^lt^  shear  noise  source  term  ) or  u,^^  ^Ux^/(3t^-v 

self  noise  source  term)  were  measured  at  various  source  positions  in  the  jet. 
Source  distributions  were  inferred  therefrom  over  slices  of  jet  normal  to  the 
jet  axis.  They  were  found  to  be  strongly  pear-shaped,  rather  than  axisymmetric , 
the  small  end  of  the  'pear'  pointing  toward  the  observer.  (This  is  not,  however 
incompatible  with  the  axisymmetry  of  far  field  sound  intensity.)  Self  and  sheai 
noise  spectra  have  been  constructed  from  the  measured  cross-spectral  densities 
by  a method  consistent  with  the  postulated  self/shear  noise  formalism.  The  two 
spectra  exhibit  comparable  amplitudes  and  virtually  identical  shapes,  but  are 
displaced  substantially  in  frequency:  all  this  is  predicted  by  the  theory. 

Self  and  shear  noise  spectra  extracted  from  far  field  jet  noise  intensities  via 
an  algorithm  of  Nossier  and  Ribner  exhibit  the  same  behavior.  On  the  whole 
both  sets  of  spectra,  although  derived  from  vastly  different  experimental 
procedures,  are  compatible. 
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